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A Neutron Absorptiometer for the 
On-Stream Analysis of Boron 


CARROLL L. PLEASANCE 
Process Control Development Engineer* 


INTRODUCTORY 


Carroll L. Pleasance graduated from the University of Washing- 
ton in 1951 with a B.S. Degree in Electrical Engineering. Since 
graduation he has been employed by the General Electric Company 
at its Hanford Atomic Products Operation. Since 1954 he has 
worked on the development and application of control instrumenta- 
tion for chemical process plants. While working in Richland, he 
has been enrolled at the University of Washington Center for 
Graduate Study at Hanford, where he is working toward the degree 
of Master of Science in Engineering with a Nuclear Option. 

The paper presented here marks an occasion because Mr. Pleasance 
is a candidate for one of the first degrees to be obtained through 
academic work mostly done at Hanford and because it is the first 
work by a Hanford student to appear in Trend. Mr. Pleasance’s 
paper deals with research work he performed as part of his assigned 
professional duties at the Hanford project; he intends to submit a more detailed report of this experimental work in 
fulfillment of the thesis requirement for the M.S. degree that he expects to receive this spring—KERMIT B. BENGTSON, 


C. L. Pleasance K. B. Bengtson 


Director, University of Washington Center for Graduate Study at Hanford. 


Hanford Atomic Products Operation has under- 
taken development studies for reprocessing power- 
reactor fuels. This program involves the chemical 
extraction and recovery of valuable nuclear materials 
from the fuel of privately owned power reactors. 
These fuels contain higher concentrations of Uran- 
ium-235 than is routinely processed in the existing 
separations plant ; and the presence of this Uranium- 
235 in solution may create high neutron fluxes in 
process vessels. The possibility exists of producing 
critical mass in vessels that are not critically safe. 
Solution concentration control is used as a primary 
means of criticality control. However, boron is to be 
added to some process solutions to increase the ther- 
mal neutron cross section: this will reduce the neu- 
tron flux and avoid the possibility of an uncontrolled 
critical reaction. Thus a boron addition will afford an 
extra margin of nuclear safety during processing of 
power-reactor fuels. 

An on-stream analyzer is needed to determine 
whether process solutions contain sufficient boron to 
assure the extra safety factor. This paper describes 
such an analyzer, an absorptiometer. It was designed 
for installation at a particular point in the separa- 


*The author works for Hanford Laboratories, General 
Electric Company, Richland, Washington, where the work 
described in this paper was performed under contract No. 
AT-(45-1)-1350 for the U.S. Atomic Energy Commission. 
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tions plant. Shown in Fig. 1 is the physical arrange- 
ment of the operating gallery, pipe gallery, and pro- 
cess cell. In the operating gallery, cold chemical solu- 
tions are prepared in chemical make-up tanks prior 
to addition to the process. These solutions flow by 
gravity from the operating gallery through the pipe 
gallery to the process cell. The analyzer will be in- 
stalled in the pipe gallery; readout will be done in 
the operating gallery. 

The solutions that pass through the analyzer are 
60% nitric acid, 50% sodium hydroxide, 38% sodium 
nitrate, and water. Each of these solutions contains 
approximately | g per liter of boron in the form of 
boric acid. Owing to the large volumes of solution 
that must be handled, the analyzer should allow a 
flow of at least 50 gal per minute to the process with 
a full make-up tank. The instrument apparently did 
not need to be more precise than +10%, the idea 
being that when it detected too low a boron concentra- 
tion the boron would probably be completely absent 
from the solution. The instrument must respond and 
indicate a lack of boron in solution in less than 15 
sec. This installation is to be a prototype for plant 
evaluation; therefore only alarm circuitry, such as 
indicating lights, is necessary. Control circuitry may 
be added when the monitor has demonstrated satis- 
factory performance. 

Stromatt surveyed the literature on methods of 
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analysis for boron. His study reveals that boron 
may be analyzed by colorimetric, volumetric, spectro- 
graphic, flame photometric, polarographic, polari- 
metric, and neutron-absorption methods. Most of 
these methods are not readily adaptable to rapid on- 
stream analysis of boron since they generally measure 
total boron rather than selectively analyze for the 
neutron-absorbing isotope boron-10. The thermal 
neutron cross section of the solution is the specific 
property that needs to be determined. Since neutron- 
absorption techniques measure precisely this property, 
the development of a neutron absorptiometer was 
begun. 


The Sensing Unit 

A coaxial arrangement of the sensing unit was 
planned as shown in Fig. 2. DeFord and Braman 
found that such a configuration produces the best 
counting geometry.” The sensing unit consists of 
three parts—a Source and Moderator, a Sensing Cell, 
and a Seal Pot. 


Source and Moderator. The source used in the 
monitor is a plutonium-beryllium neutron source 
emitting 1.5 x 10° neutrons per sec. The source pro- 
duces neutrons by an alpha — neutron (a, m) reac- 
tion, the alpha particle being produced from the plu- 
tonium decay: 


+ — + on’ + Q. 


The energy associated with the neutron is usually 
greater than 5 mev; so this is called a “fast” neu- 
tron. To be utilized for measurement of boron con- 
centration, the neutrons must be slowed to “thermal” 
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energies (less than 0.05 ev) to permit capture by the 
boron in the solution undergoing analysis. This 
slowing down is accomplished by placing a paraffin 
moderator around the source. The fast neutrons 
enter the moderator and undergo elastic collisions 
with hydrogen and carbon atoms. This both scatters 
the neutrons throughout the moderator and slows 
them to thermal energy levels. Thus there is a three- 
dimensional distribution of thermal neutrons traveling 
throughout the moderator. A certain percentage of 
these neutrons travel through the sensing cell and are 
measured by the neutron detector tube. As the 
thermal neutron cross section of the solution in the 
sensing cell increases, the counting rate of the 
detector decreases; therefore the counting rate is 
inversely proportional to the thermal neutron cross 
section of the solution. 

Sensing Cell. Once the counting geometry of a 
particular configuration of source, moderator, and 
sensing cell has been determined, the calibration 
curve may be approximated by calculating the change 
in macroscopic cross section as the solution concen- 
tration in the sensing cell is varied. For a high sensi- 
tivity to low concentrations of boron, the sensing cell 
and moderator housing must be fabricated in such 
a manner and of such materials as to minimize the 
absorption cross section. Experimentally it was de- 
termined that a thin-walled stainless-steel cell with a 

(Continued on page 31) 
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Water Disinfection at a Natural Swimming Beach 


ROBERT O. SYLVESTER 
Professor of Sanitary Engineering 


In February, 1959, the 
University of Washington, 
through its Department of 
Civil Engineering and at 
the request of the City of 
Seattle, entered into a con- 
tract with the Seattle Park 
Commissioners for a study 
of Green Lake that would 
result in recommendations 
for rehabilitation of the 
Lake. The Lake has a sur- 
face area of 256 acres, a 
shoreline of 2.9 miles, and 
an average depth of 12.5 feet. It contains some 1088 
million gallons of water. It is surrounded by a 
thickly populated residential area in the north-central 
portion of Seattle. Green Lake lies on the margin of 
the University District. 

Because of its location and character, Green Lake 
has a present and potential recreational usage that 
is not excelled by any other lake of its size in the 
United States (see Figs. 1 and 2). The full realiza- 
tion of this potential has been jeopardized by heavy 
algal blooms during the swimming season, by “swim- 
mers’ itch,” sewage backflow into the lake, and by 
the maturation of the Lake to the stage of rapid 
shoaling, with the appearance of swampy character- 
istics about its shoreline. During the period of 1953- 
57, the yearly average number of swimmers was 
about 350,000. Green Lake was closed to swimming 
during a portion of the 1958 season because rela- 
tively high coliform bacterial densities were observed 
in the Lake water and because the North Trunk 
Sewer had previously collapsed and backflowed into 
it. 


R. O. Sylvester 


A portion of the contractual agreement for the 
overall Lake study called for an immediate study of 
what interim measures might be taken to improve 
water quality at the two Green Lake bathing beaches 
in time for the 1959 season. This article summarizes 
the report on these interim measures,’ including the 
operating data for the 1959 season. The report on 
the general lake study, with recommendations for a 
program of lake rehabilitation, was submitted to the 
City of Seattle in February, 1960.? 


Physical Characteristics of Beach Areas 


Of the two swimming beaches on Green Lake, one 
on the west shore and the other on the east, the West 
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Beach has the best facilities and the most attractive 
surroundings (Fig. 3). Each beach has a supervised 
frontage of 380 ft and an offshore limit of 145 to 200 
ft, giving each a water area of about 72,000 sq ft. The 
water at the West Beach extends to a depth of 14 ft, 
and the water volume in the beach area is 5.2 
million gal; the East Beach (Fig. 4) extends to 20 
ft and has a water volume of 6.3 million gal. Most 
of the West Beach has a sand and gravel bottom 
except for the very outer portion, where mud depths 
range from 3 to 13 ft. At the East Beach, the entire 
bottom is mud except for the central portion and the 
shallow-water areas, which are sand. Mud depths 
range from 7 to 13.5 ft throughout the swimming 
portion of the beach area. 

Water currents were measured at both beaches by 
using drogues and dye. Light winds produced a 
water movement of 150 to 220 ft per hr in the upper 
6 ft, and moderate winds produced velocities of 230 
to 300 ft per hr. On a calm morning, several hours 
after a moderate wind had been blowing, erratic 
water movement still persisted with velocities of 90 
to 125 ft per hr. The average surface wind velocities 
and direction at the Seattle Weather Station during 
the swimming season (for 20-yr and 38-yr averages) 
are as follows: 


June July August September 


Avg wind velocity, mph 42. 39 Ks 3.9 
Predominant direction Ss N N Ny 


Thus in summer, with light northerly or southerly 
winds, we could expect that the water in the swim- 
ming areas would be exchanged by surrounding lake 
water in from two to four hours. During periods of 
calm weather (with no appreciable wind for several 
days) the water in the bathing area would experience 
only a slight exchange with the surrounding water. 
These periods of prolonged calm are infrequent. 


Schemes of Water-Quality Improvement 


Any water-treatment measures that might be used 
to improve the water quality for swimming in Green 
Lake during the 1959 season needed to have public 
acceptance of their reliability at the outset; be 
capable of simple operation; be limited in their 
effect on the Lake; be relatively inexpensive, since 
they might not fit into a program of long-range lake 
improvement; be suitable for easy construction, to 
permit their completion in time for the swimming 
season ; and, above all, needed to eliminate the health 
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hazard that has been associated with swimming in 
Green Lake. 

A number of schemes that might have resulted in 
improved water quality at the swimming beaches 
were investigated, including water treatment with a 
recirculation system of lake water, the addition of 
copper sulfate for algae control to improve the 
water’s appearance and the functioning of a purifica- 
tion system, diking off the swimming areas, admission 
of city water to the beaches, and chlorinating the lake 
water at the beaches. 

The installation of a system at the bathing beaches 
to recirculate lake water through a chemical precipi- 
tation plant, with chlorination to follow, was not 
practicable because of its high cost and because the 
water opposite a beach is being continually exchanged 
with outer lake water, so that treatment of the 
entire lake becomes necessary. A water-treatment 
plant capable of treating the average water supply 
for a city of 20,000 persons would be required just to 
exchange all the lake water through the plant but 
once a year. Such a plant would also have no 
appreciable effect on the overall algal problem of the 
lake unless large quantities of lime and alum were 
added to precipitate phosphorus. 

An alternative remedy would be to construct a 
dike around each of the two bathing beaches, dredge 
the bottom sediments, and replace them with sand. 
A dike constructed of a nonorganic, silt-free mixture 
of sand and gravel would require frequent repair 
unless protected with concrete or similar walls. 
Recirculation of this contained water through a filtra- 
tion plant would be required, since dilution and 
dispersion of this diked-in water with outside lake 
water would now be precluded. If these diked-in 
areas were to contain as much water as is now 
available at the beaches (five to six million gal each), 
and if State Health Department requirements for 
water exchange and filtration rates were observed 
(exchange every 6 hr with a filtration rate of 3 gpm 
per sq ft of sand in a pressure filter), quite an array 
of equipment would be required at each beach: two 
7,000-gpm recirculation pumps; one alum and one 
lime feeder, each capable of feeding 2500 Ib per 24 
hr; one chlorinator with a feeding rate of 250 lb per 
24 hr; and eighteen 10- by 25-ft pressure filters. The 
expense of constructing and operating this system 
to accommodate the same number of bathers now 
served at the natural beaches was not feasible. If the 
water volumes were reduced at these beaches to about 
half their present volumes, the cost would still be 
very high; and the remainder of Green Lake would 
be no better than before for the many other recrea- 
tional uses. 

The best scheme for meeting the requirements 
appeared to be a system that would diffuse chlorin- 
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ated water into the swimming areas rapidly enough to 
overcome the chlorine demand of new lake water 
constantly being brought into the beach area by wind- 
driven currents. Treatment at the West Green Lake 
beach was recommended over the one on the east 
shore because the beach bottom, water depths, and 
water currents are more suitable for proper function- 
ing of a disinfection system and because the East 
Green Lake area has, besides the beach, a swimming 
pool, baseball field, and boat-rental facilities. Another 
recommendation was that the East Green Lake beach 
be opened for swimming during the 1959 season, 
to be closed only if there was evidence of undesirable 
fecal pollution from human sources. Much valuable 
data on the effectiveness of the treatment proposed 
and on the epidemiological relationship* between the 
two beaches could be obtained if both beaches were 
kept open. 

Disinfection of natural bathing areas was known 
to have been used at the following locations, most of 
which are ponds having little water motion: 


Blue Lake near Portland, Ore. 

Wallace Lake near Cleveland, Ohio. 

Watchung Lake, Plainfield, N.J. 

Skytop Lodge, Barret Township, Pa. 

Sons of Poland Camp, Inc., Morris County, N.J. 
Field’s Grove Pool, Nashua, N.H. 

Community Park Pool, Woodbury, N.Y. 
Fountain Spring Lake, Wanaque, N.J. 
Riverdale Pool, Riverdale, N.J. 

Paramus Beach, Paramus, N.J. 
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Required Chlorine Dosages 

The quantities of chlorine required to kill bacteria 
and to maintain a chlorine residual of at least 0.2 ppm 
in the lake water were determined by pilot-plant 
experiments with water diffusers at the Aqua The- 
ater Pool (another Green Lake recreational attrac- 
tion) and by laboratory experiment, in which lake 
water was “doctored” in an attempt to simulate the 
extreme conditions of water quality that might be 
found at a bathing beach with a crowded bathing 
schedule and with a heavy sediment and algal load. 
(Data were not available on water-quality constitu- 
ents present at a Green Lake bathing beach during 
a summer swimming season). The results of these 
laboratory and field tests, together with those ob- 
tained from the Seattle Water Department, are given 
below. These values gave a rate of chlorine feed 
ranging from 420 to 1670 lb per 24 hr when applied 
to the West Beach requirements. Assuming a 9-hr 
day of chlorination, the actual amount of chlorine 
added to the water at the West Green Lake bathing 
area would range from 160 to 320 lb per day. The 
following tabulation shows the dosage required to 
give a chlorine residual of 0.2 ppm within the stated 
periods of time: 


Detention 0.5 2 3 4 5 
Times (hr) 
o o o 
Chlorine 1.4/0. 
Req. (ppm) 2.4 4.9 5.3 5.6 6.9 


Details of Disinfection System 


The West Green Lake disinfection system was con- 
structed to function as follows: 


1. Three headers containing orifices, perpendicular to the 
shore, inject chlorinated water into the swimming area 
(Figs. 2, 5, and 6). When a northerly wind is blowing, 
producing a southerly current, the south header is closed 
off. When a southerly wind is blowing, the north header 
is closed off. Under conditions of no wind (little current) 
all headers are operated. The quantity of chlorinated 
water injected into the different areas is controlled by 
varying the chlorine feed rate and by use of the valves 
shown in Fig. 6 (pump-house plan). The amount of chlo- 
rine to be fed and the proper setting of the control valves 
is determined by measurements of chlorine residual in the 
water, the desired residuals having been determined from 
bacteriological tests. 


2. Chlorine is distributed by the recirculation of lake water 
through a pump at the rate of about 500 gpm, the quantity 
depending upon the number of headers in use. City water 
would have been more desirable for this use as it has little 
chlorine demand, has adequate pressure to eliminate the 
need for a pump, does not tend to clog the orifices, and 


APRIL, 1960 


PLES AT 1S’ CENTERS 


ORIFICES ALTERNATING FROM HOR TO BELOW HOR -7" cenrens 
2 ORIFICES MORIZONTAL ORIFICES FACE 


ELEV «15450 


SOUTH HEADER 


ba 30° NORTH HEADER:- 
174” ORWLED HOLES 
CENTER HEADER: - 
SIMILAR EXCEPT 22 ORIFICES FACE BOTH DIRECTIONS IN 4” PLASTIC PIPE, WHICH 
1S 6 BELOW WATER SURFACE. FIRST 4 ORIFICES ARE HOR, 2ND 4 ALTERNATE 
FROM 16° BELOW TO HOR. ALL OTHERS ARE 30° ABOVE HOR. 


SIMILAR EXCEPT ORIFICES FACE SOUTHWARD AND ARE 
SPACED 6 CENTERS. FIRST 4 ORIFICES ARE HOR. 


TROLLEY BM OVER © 
| 
ar 
~ 
2 at H 
10 Cr ORATOR) 
PLASTIC CHLORINE SOL LINE 4) 18 mec 
PLAN-PUMP HOUSE 


SCALE-3/:6" FT 10 PUMP HSE 
PLAN- CHLORINATION. 


Fic. 6. DETAILS OF DISINFECTION LAYOUT 


would give some desirable dilution to the lake water. 
However, lake water was recommended because of the 
possible shortage of city water if the summer was a dry 
one and because the cost of buying, installing, and oper- 
ating a recirculation pump would be less each year than 
using city water (considering the commercial value of the 
city water). When the pump is operating efficiently, with 
the discharge valve or valves wide open, the discharge and 
gage pressure should be about as follows: 


Headers Working Gage Press.( psi) Pump Disch.(gpm) 
South or North Only 80-84 380 
Central Only 65 530 
South or North 

and Central 55-60 580 
South and North 65-70 510 
All Headers 40 720 


3. The chlorinator is a manual machine having an operating 


range of 200 to 2000 Ib of chlorine per 24 hr. City water 
at the rate of about 50 gpm is used to operate it and to 
inject the chlorine solution into the pump suction through 
2-in. plastic pipe. Chlorine is supplied from 6 one-ton con- 
tainers arranged in parallel on a central header. The use 
of these large containers in parallel made the purchase of 
an evaporator unnecessary. (A maximum of 400 Ib of chlo- 
rine per day can be withdrawn from a ton cylinder without 
the use of an evaporator, which costs about $2300.) An 
empty cylinder alarm was used instead of the costly scales. 
Arrangements were made with the City Water Depart- 
ment to deliver the chlorine cylinders since they are 
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SAMPLING STATIONS 


identical to those obtained by the Water Department on 
contract. The chlorinator was started each morning in 
sufficient time for the build-up of 0.3 ppm of chlorine 
residual before swimming commenced. A separate build- 
ing was erected for the chlorinator and chlorine cylinders 
as indicated in Figs. 5 and 6. 


4. A 500-gpm recirculation pump (an iron-body turbine type 
with 160-ft head) powered by a 220/440-v electric motor 
was installed in a pump and control house located in front 
of the bathhouse (Fig. 5). The pumphouse floor was de- 
pressed 5 ft below the ground surface to reduce the build- 
ing’s visible height and the length of pump shaft required. 
Electric conduits to the pumphouse were buried. Three 
plastic valves to control the chlorine dosage to each distri- 
bution line are located in the pumphouse. A 6-in. pump 
suction line extends into 10 ft of water between the diving 
platforms. 


5. The chlorine diffuser lines or headers are plastic pipes 
supported from piling driven on 15-ft centers. This type 
of support permits easy removal of the pipe for cleaning, 
adjustment, or repairs. The headers on the north and 
south boundary of the bathing area are located 30 in. 
below the water surface; the center header is 6 ft below 
the water surface where it is opposite the diving plat- 
forms. Orifice details for the headers are shown in Fig. 6. 
To permit flushing for cleaning, each header has a valve 
at its outboard end. All pipe on land and near the water’s 
edge is buried. Signs are posted on the piling asking swim- 
mers to stay off the pipe and to keep away from the piling. 

6. The complete equipment and construction cost for the fa- 
cility was $33,949. 


Operating Requirements 

The chlorination facility was placed in operation 
on June 16, 1959, with control of the operation under 
the head lifeguard and the Evan’s Pool maintenance 
personnel. The arrangement did not prove satisfac- 
tory; so a full-time operator (a University student) 
was hired on June 24 for the remainder of the swim- 
ming season. This operator had these duties: (1) 
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starting the pump and chlorinator about 1 hr before 
the beach opened so that a chlorine residual of ap- 
proximately 0.3 ppm could be built up prior to the 
commencement of bathing and swimming; (2) col- 
lecting chlorine-residual samples every 42 to 2 hr, 
depending upon the conditions, at 6-in. and 6-ft 
depths from the 7 stations shown in Fig. 7; (3) ad- 
justing the chlorinator feed rate and header-control 
valves to give the desired chlorine residual in the 
beach areas where there were swimmers and bathers ; 
(4) collecting periodic samples for bacteriological 
analysis; (5) keeping a running record of the opera- 
tion as shown in Table I; (6) maintaining the sys- 
tem by flushing the headers, pumping out the pump 
pit, and changing chlorine cylinders; and (7) shut- 
ting the system down at the close of the swimming 
day. 

The total cost of about $1930 for the 1959 sea- 
son included operation, $1100; chlorine, $650; and 
power, $178. 


Operating Conditions 

The amount of chlorine required for disinfection 
varies with the sections of the beach area being used, 
the wind direction and velocity, the amount of sun- 
light, the water quality, and the number of bathers 
in the water. The lifeguards limit the bathers and 
swimmers to designated areas, depending upon the 
number present. With 20 or fewer bathers, they 
are limited to the south end of the south section; 
with 20 to 100, they are limited to the south section ; 
and with over 100, the entire beach area is opened. 
Swimming classes are held inside the rope along the 
north section. Only the south header need be oper- 
ated when the wind is southerly and just the south 
section is in use. If there is no wind, the central or 
the central and south headers should be operated. 
A northerly wind usually means a good day for 
swimming, and the entire beach area is open; then 
all three headers or the north and central headers 
are used. 

Table I is a sample from the record sheets that 
summarize the chlorination plant operation and the 
chlorine residuals obtained in 1959. This sampling 
from the summary sheets for 72 days shows the 
record for 14 days during the peak of the bathing 
season. On the average, about 60 separate chlorine- 
residual readings were taken each day. After discov- 
ery that the central header had been installed upside 
down, it was realigned on July 1. The north header 
outboard of the fifth piling (last half) is turned so 
that the orifices all point downward instead of alter- 
nating between horizontal and 15° below horizontal. 
This misalignment should be corrected before the 
next swimming season. With the realignment a gate 
valve should be installed about midway in this header 
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TABLE I 


SAMPLE FROM WEST GREEN LAKE 
CHLORINATION SUMMARY (JULY 26-AUG. 8, 1959) 


mg/l1 
No. of Cl, Feed Ch Opera- Avg. Chlorine Residual at Sampling Sta.t Headers in 
Date Wind | Bathers* Range Used | tion (hr.) Operation 
(Lb.) (Lb.) 1 2 3 4 5 6 7 
7-26-59 | NW 20-300 | 4-500 2057 9 0.45 | 0.38 | 0.2 10.6 10.7 10.8 104 Pe 
27 Rain-S | 0-20 2-600 2202 9 0.2 0.15 | 0.1 0.2 0.35 | 0.2 0.1 S 
28 NW 20-250 5-900 2457 9 0.1 0.1 0.1 0.5 0.7 0.4 | 0.3 | 
29 NW 10-200 | 5-600 2663 9 6:2 10.15 1 
30 NW-N 50-400 5-900 2949 8.5 0.2 0.1 0.1 0.75 |}0.6 | 0.4 | 0.2 N&C 
31 W-NW | 60-500 7-1000 3259 9 0.3 0.2 0.2 1.0 | 0.4 | 0.25 | 0.15 | NC&CS 
8-1-59 | S 0-150 2-700 3443 9 0.15 | 0.1 0.1 0.1 0.75 | 0.4 | 0.2 CS&S 
2 S-SW 0-50 2-700 3651 8.5 0.2 ©.15 | 4.15 1¢@.5 0.75 | 0.25 | 0.15 | CS&S 
SW 0-150 3-800 3911 9 0.2 0.1 0.1 0.15 10.7 0.3 0.2 Ss 
- NW 0-10 2-700 4061 9 0.05 | 0.05 | 0.05 | 0.2 0.5 0.2 0.15 | C 
5 Ww 20-200 | 6-700 4311 9 0.5 | 0.35 | 0.3 10.8 10.5 1.03 
6 NW-N 20-300 6.5-1100 4607 9 0.2 0.25 | 0.2 0.7 0.75 | 0.5 0.3 NC&C 
7 NW-N 30-200 6.5-1100 4898 9 0.8 |0.4 |0.4 | 0.7 0.75 | 0.5 0.2 NCS&C 
8 NW-N 5-400 7-800 5167 9 0.4 | 0.3 0.3 0.6 | 0.75 |0.6 | 0.3 Cc 


* Number of bathers shown is the actual count of those in the water (range) at any one time during the day. 
+ The number of chlorine residuals taken daily from each sampling station varied between 6 and 16. Stations 2, 3, 6, and 7 were sampled from 


near the surface and at mid-depth. Chlorine residual was 


lowed to decline at stations 1, 2, and 3 when the lifeguards restricted swimming 


to the south portion of. the beach. This was the case during periods of low swimming population. 


so that chlorine can be shut off in its outer half when 
there is no swimming beyond the ropes in this sec- 
tion of the beach. 


A total of about 15,200 lb of chlorine was used 


in 1959 for 72 days of swimming, or an average of 
211 lb per day. The rate of chlorine feed ranged 
from 200 to 1100 Ib per 24 hr; an average operating 
day was about 8.5 hr. 

On one occasion the intake screen was lifted off 
the end of the suction pipe by a group of boys; and 
shortly thereafter the pump lost its suction, chlorine 
continued to feed into the system, and the odor of 
chlorine became pronounced in parts of the swim- 
ming area. After shutting down the system, police 
skin divers replaced the screen, and the facility was 
returned to operation. The recirculation pump later 
developed a serious vibration, and it was shut down 
on September 2. After its removal, the pump im- 


peller was found to be seriously corroded by the 
chlorine fed into the pump suction when the pump 
was not pumping water. An automatic switch and 
alarm to be installed will shut off the pump and the 
chlorinator if ever again the pump loses its suction. 


Operating Results 


The chlorination facility operated very well with 
the exception of the troubles mentioned in the pre- 
ceding section. There were no complaints of eye 
smarting or chlorine odors. A few dozen trout were 
killed on the first day of operation, but no kills of 
any significance were reported thereafter. The West 
Green Lake Beach was one of the most popular in 
King County. As shown in Table II, it had the best 
bacteriological record of the six beaches in King 
County sampled during the summer of 1959 by the 
King County Health Department. Coliform bacteria 


TABLE II 


BATHING-BEACH BACTERIOLOGICAL DATA* 
KING COUNTY, JUNE 25 - SEPT. 2, 1959 


GREEN LAKE WASHINGTON Pucet Sounp LAKE 
GOLDEN SAMMAMISH 
West Beacu | East Beacu | SEWARD PARK JUANITA GARDENS STATE PARK 
No. of values 113 105 113 113 112 113 
Median valuet 230 230 “60 230 60 
Range of values 0 to 2400 0 to 24,000 | 0 to 13,000 0 to 7000 0 to 7000 0 to 2400 
Avg No. bathers 130 104 110 110 43 114 
Avg water temp. 70.3 71.1 69.0 70.5 58.2 71.9 
Coliform per bathert 0.18 4:3 a4 0.54 5.4 0.53 


* Most probable number of coliform bacteria per 100-ml sample, computed from Seattle-King Co. Health Dept. data. 


t Less than 45 proportioned equally among 0, 15, 30, and 45. 


Median 


t Coliform per bathers 
Avg No. Bathers 
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TABLE III 
STATISTICAL ANALYSIS OF GREEN LAKE 
BACTERIOLOGICAL DATA, JUNE 25-SEPT. 2, 1959 


ARITHMETIC GEOMETRIC STANDARD CoeEF OF VARIATION | RANGE CONTAINING 

SouRCE MEAN MEAN MEDIAN DEVIATION or Geo. MEAN 50% or VALUES 
West Beach* 1054 30 23 4.48 0.15 11-83 
West Beacht 20.7 2.21 2 11.35 5.14 0-13.6 
West Beacht 7.06 3.94 4 4.11 1.04 1.48-10.2 
East Beach* 10114 196 230 6.74 0.32 55-720 
East Beacht 27.9 16.2 14 3.72 0.23 11. 25-23.4 
East Beacht 9.81 4.06 3.6 12.0 2.96 0-21.8 
Green Lake§ 17.5 3.16 3.6 14.41 4.56 0-19.1 


Kin 


U. of Winkincon data: 67 samples at W. Beach, 38 samples at E. 
} U. of Washington data: 30-39 samples obtained just outside bathing area and up wind. 
§ U. of Washington data: 94 samples obtained at regular sampling stations—all Calle bathing areas. 


{ Values influenced by a few extreme coliform densities. 


concentrations were very low, the median being less 
than one-tenth of the tolerable limit of 240* estab- 
lished by the State Health Department for Wash- 
ington water. A normal tolerable limit of 1000 is 
used in California and many other states.** If the 
average chlorine residual can be kept between 0.4 
and 0.6 ppm during the next swimming season, the 
coliform bacterial density in the bathing and swim- 
ming areas should be even better. It should be noted 
that the coliform bacteria are merely indicators of 
water quality. 

The following comparison, which applies to Table 
II, points up the record for the six King County 
beaches mentioned in the preceding paragraph: 

1. The West Green Lake Beach had the highest water qual- 
ity (bacteriological) and the most bathers. 

2. The East Green Lake Beach (unchlorinated) had coli- 
form concentrations greater than Juanita Beach on Lake 
Washington and greater than the Lake Sammamish State 
Park. It was comparable with Seward Park on Lake 
Washington and Golden Gardens on Puget Sound. 

3. Golden Gardens bacteriological data are influenced by the 
nearby North Trunk Sewer Outfall and the discharge 
from the Ship Canal. 

4. The water temperature at all beaches was about the same 
except for Golden Gardens, where it was 12° cooler. 
Table III presents a statistical analysis of coli- 

form bacteria densities observed at the two Green 

Lake bathing beaches by the University and by the 

Seattle-King County Health Department during the 

1959 bathing season. Health Department samples 

were collected in the area of greatest swimming con- 

centration, three samples being taken (surface, mid- 

dle, and bottom) at knee depth and a similar three 

at a depth of 5 ft. University samples were collected 
from 6-in. and 6-ft depths in three locations—outside 
and upwind of the swimming area (control samples), 


* These limiting values usually apply to bacteria from 
humans rather than from wild life. But the bacteria found 
in Green Lake during the 1959 season were largely from 
wild life; and they have a different sanitary significance 
from those that come from humans. 
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H. D. data: 113 samples at W. Beach, 105 samples at E. Beach. 
Beach. 


in the wading area, and around the platform where 
most swimmers were congregated. Table III also 
shows the bacteria densities obtained from through- 
out the Lake (exclusive of bathing areas) for the 
same time period. The data in Table III are the 
basis for a comparison of the chlorinated water at 

West Beach and the untreated water at East Beach: 

1. The bacteria concentration was sufficiently low for swim- 
ming at both beaches. 

2. The bacterial densities at the West Beach were one- 
seventh to one-tenth of those at the East Beach. 

3. The water quality inside the swimming area at the West 
Beach was about the same as that outside the swimming 
area, despite the presence of many bathers. 

4. The coliform concentration at the East Beach was about 
three to four times as high inside the swimming area as 
outside it, though the count for the inside area is not high. 
The bathers contributed these additional bacteria. 

5. The water quality in the swimming areas did not differ 
markedly (according to University data) from that of the 
lake water in general, indicating the large water exchange 
brought about by wind-induced currents. 

6. University data gave lower values than did the Health 
Department data because of a difference in the sampling 
points. 

7. The statistical evaluation was hampered by a large num- 
ber of zero values and by many Health Department values 
that were recorded as “less than 45.” 


Recommendations Made to the City 
Now that the North Trunk Sewer has been re- 
paired, it is unlikely that any more sewage will be 
discharged into Green Lake. The Lake has also puri- 
fied itself of former sewage discharges. However, a 
great many persons have become skeptical of the 
water quality of the local beaches and are fearful 
about swimming in any natural body of water. Oth- 
ers will, of course, swim anywhere at any time. By 
providing one natural bathing area with a controlled 
bacteriological quality, the City would be making a 
truly worthwhile contribution to the general well- 
being of these many people who wish to swim in a 
natural setting but are hesitant to do so because of 

(Continued on page 32) 
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On the Brittle Fracture of Structures 


PAUuL C. Paris 
Assistant Professor of Civil Engineering 


Galileo in 1638 gave a 
lucid description of the 
character of internal forces 
that resist separation with- 
in material bodies.t Since 
his time only two signifi- 
cant types of separation or 
fracture theories have been 
proposed. The first is the 
concept that the “state of 
stress” at a point may be 
regarded as the cause of in- 
cipient fracturing. Coulomb 
in 1776 made use of this 
idea to describe material failures. This concept has 
been carried down to date as a criterion of yielding 
at a point in a body; however, it has not been suc- 
cessful in describing the “brittle failures” of many 
structures over the past century. 

The second theory proposes that all materials and 
structures contain flaws and that under load the 
worst of these flaws will grow, causing fracture. 
Griffith in 1920 developed the first successful theory 
of flaw growth in brittle solids.? His methods have 
been modified and developed as a criterion of “brittle 
fracture” of ductile materials, especially by Irwin.* 
These relatively new concepts are presently em- 
ployed by the aircraft industry and others in ana- 
lyzing fractures in structures composed of ductile 
metals. The choice of materials for aircraft struc- 
tures is based in part on the results of these modern 
fracture analyses. Hence this article will attempt to 
describe these theories and to illustrate a few appli- 
cations, as well as to indicate some of the current 
trends in fracture analysis. 


Griffith Theory Applied to Ductile Materials 


Griffith derived a stability criterion for flaws 
(cracks) in ideally brittle materials such as glass. 
He reasoned that a crack in a plate of material will 
grow if the energy supplied from the stress field, @ , 
exceeds the surface energy of the new crack surface, 
2T,* during an increment of. extension. Hence for 
the commencement. of instability, 


G=2T. (1) 


Using stress analysis to compute G for the config- 


P. C. Paris 


* A 2 is required with T since the surface area of both 
faces of a crack must be considered. 
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WILLIAM E. ANDERSON 
Research Engineer, The Boeing Airplane Company 


uration shown in Fig. 1, he 
found that 


or 
-(3) 


Now, the right side of 
Eq. (3) is a material con- 
stant; hence flaw stability 
depends on both stress level 
and flaw size. 

The tip of a crack has a theoretical stress-concen- 
tration factor approaching infinity; hence in: ductile 
materials a zone of yielded material always accom- 
panies a crack tip. As the crack extends, plastic de- 
formation occurs that dissipates energy. This adds to 
the surface energy, 27. 

Equations (1) through (3) can be modified by 
the additional term, dP/dA, the plastic work per unit 
of crack surface area created in crack extension. From 


W. E. Anderson 


Fic. 1. THE GRIFFITH CONFIGURATION 
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FiG. 2. TYPICAL CRACKED-STRENGTH BEHAVIOR 
these considerations, Eq. (3) becomes 


Experiments and additional theoretical considera- 
tions show that dP/dA is a constant for a given thick- 
ness of a material provided that the nominal stress 
level, o, is below the yield point, o,,. Moreover, 2T 
is negligibly small compared to dP/dA in ductile 
metals. 

Figure 2 shows Eq. (4) plotted for a typical ma- 
terial. By using severe yielding as a criterion of 
failure also, the curve is cut off at the yield point of 
the material. If a structure contains flaws larger than 
the transition size, brittle fracture occurs rather than 
gross yielding when the structure is overloaded. Such 
behavior is undesirable in that the structure should 
be able to dissipate the energy of occasional over- 
loads without the occurrence of fracturing. More- 
over, the Theory of Plastic Design indicates that the 
ability of structural components to yield rather than 
to fracture allows a redistribution of forces in redun- 
dant structures that enhances their strength. 


So-Called High-Strength Metal Alloys 


The common alloys of aluminum used in aircraft 
structures are numbered 2024-T3, 7075-T6, and 
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Fic. 3. COMPARISON OF COMMON ALUMINUM ALLOYS 


7178-T6. The latter two are often called “high- 
strength” alloys because their tensile yield strengths 
are half again higher than 2024-T3. However, their 
dP/dA values are considerably less than 2024-T3. 
Curves of their behaviors that have been prepared 
(in the manner of Fig. 2) are superimposed in Fig. 
3. (A sheet thickness of 0.064 was selected as a 
typical gage for Fig. 3, and the required data were 
furnished by the Boeing Airplane Company.) 

As can be seen in Fig. 3, the so-called high- 
strength alloys are stronger only if flaws are very 
small. Since accidental damage or fatigue cracks 
are likely to occur during the service life, 2024-T3 is 
a much safer alloy for airframes. 

Thus modern trends toward higher-strength metal 
alloys in all types of structures will require special 
considerations of the possibility of brittle fracture. 
Special tests must be performed to evaluate the tear 
resistance (i.e., dP/dA) of materials for such judg- 
ments,, since tensile stress-strain tests supply only 
the “no flaw” properties. 


The Role of the Structural Engineer 

Traditionally, the metallurgist is consulted when 
brittle failure occurs. However, since structural 
details have a great influence on brittle behavior, the 
design responsibility of the structural engineer cannot 
be neglected. For example, consider the behavior of 
a small crack emanating from the corner of a cutout. 

If the cutout is unreinforced, as shown in Fig. 4, 
the crack length, 2a, which should be used in Eq. (4) 
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for analysis, should include the width of the cutout. 
Hence very small cracks can cause brittle fracture. 
However, if the cutout is properly reinforced, the 
effective crack length can be considerably reduced. 

The tendency toward brittle fracture in a structure 

can be affected more by the design of structural 
details than any other factor. Recent ship failures, 
wherein cracks have emanated from hatch openings, 
and the failures of the deHaviland “Comet” fuselages 
by cracking from a window cutout are evidence that 
the structural engineer is only now becoming aware 
of these problems. The role the structural engineer 
must play is clear: 

1. He must obtain data from the metallurgist on the 
fracture properties (i.e. dP/dA) before making the 
choice of materials for a proposed structure. 

2. He must anticipate the possibility of flaws and cracks, 
such as fatigue cracks or damage from other sources, 


and provide structural details so that undetected flaws 
will not cause catastrophic failures. 


3. He must provide the persons servicing his structure 


with information on the flaw size or amount of 
damage the structure may safely sustain, in order that 
proper inspections may be executed. 


At present, aircraft structural design is performed 
by using this threefold procedure. The most difficult 
step is that of anticipating and analyzing cracks. 
Current research is providing stress analyses of 
many structural configurations containing cracks, 
such as was necessary to compute the left-hand side 
of Eq. (2). These analyses require the application 
of Theory of Elasticity and other advanced principles. 


re 
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Fic. 4. EFFECTIVE CRACK SIZE IN UNREINFORCED CUTOUT 
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The results of these studies are tabulated throughout 
the literature for the convenience of the design engi- 
neer, and many approximate methods of analysis are 
now available. It remains to find simple code-type 
specifications to avoid brittle fracture. 

In the future, the structural engineer will do more 
than simply obtain data from the metallurgist on 
material properties. Through a mutual understand- 
ing of the problems of design to avoid brittle fractur- 
ing, a more judicious choice of structural materials 
may be accomplished. This is to say that the 
structure itself may incorporate materials designed 
for specific applications. Hence a discussion of the 
materials-behavior aspects of the problem of brittle 
fracture is pertinent. 


Factors Affecting Tear Resistance 


The energy-dissipating capacity of a cracked 
material is the factor measured by the term dP/dA. 
This capacity is effectively due to plastic flow at and 
around the crack tip; so a study of this plasticity is 
desirable. Consideration begins on a microscopic 
level. 

In almost all materials a shear stress of some 
magnitude must occur so as to cause plastic flow. 
The operating deformation system in any given 
material usually consists of multidimensional and 
increasingly complex modes. Systems tend to resist 
further deformation all the while they are being 
deformed by the shear forces. 

Commonly used solid materials try to maintain 
contiguity with themselves; i.e., they resist the 
development of free surfaces. The deformation sys- 
tems thus tend to operate in a way that will maintain 
compatability with adjacent systems. Clearly, the 
activity of a system is influenced by the activity of 
its neighbors. 

These activities at the microscopic level are in- 
fluenced by the larger stress state imposed on them. 
Consider the region around the crack tip in a 
material of thickness involving a few tens of systems, 
all of which is in a state of plane stress. Early - 
plasticity is accomplished by initial or easy-shear 
modes. With further stretching the material thins by 
realigning of the systems. Increasingly complex 
modes can be made to operate by the shear loads im- 
posed on one system by its neighbor. The multidimen- 
sionality of the energy-dissipating systems provides 
a large amount of deformation. Associated with this 
is a large value of tear resistance. 

If the thickness of the material is reduced mark- 
edly, the number of systems is also markedly re- 
duced. Under a similar state of plane stress, the inter- 
ference of one deforming system with its neighbor is 
reduced. Fewer of the more complex deformation 
modes are achieved by the several systems. The 
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Fic. 5. OPTIMUM BEHAVIOR FOR GIVEN FLAW SIZE 


relative capacity to absorb energy is reduced and the 
thinner material will have less tear resistance than 
before. 

If the thickness is vastly increased, a state of plane 
strain occurs in the material around the end of the 
crack. This means that the material is not subjected 
to much shear. Fewer deformation systems will 
operate; and tear resistance will be lowered. 

The previous discussion is intended to be generally 
applicable to any plastically deformable solid. Knowl- 
edge of the mechanisms of deformation will allow an 
estimate of the plastic work dissipation by the given 
material in the given condition. The fracture 
mechanics of an acrylic window or a clay bank may 
not be less important than fracture of a structural 
pressure vessel, but the final discussion will be 
restricted to commercial metal alloys used in struct- 
ural applications. 


Alloy Design for Optimum Tear Resistance 


The selection of a metal for a structural function 
is normally based on cost, weight, mechanical prop- 
erties, and availability of suitable shapes. Alloy or 
processing changes can be made that will affect 
mechanical properties without significant changes in 
the other factors. This will be called alloy design. 

Important mechanical properities are yield strength 
and elastic stiffness. Gross changes in elastic stiffness 
cannot be effected, but the opposite is true for yield- 
ing. The yield strengths in some important cases can 
be altered by an order of magnitude. Since yielding 
behavior will affect tear resistances as well as other 
design considerations, a more detailed examination 
of this situation is desirable. 

Reconsider the curves shown as solid lines in Fig. 
3. If the anticipated flaw size that can be safely 
tolerated in a certain component is as shown in Fig. 
5, 7075-T6 is not safe, but 2024-T3 is too safe. 
Therefore it would be desirable to design an alloy 
that would have the property shown by the heavy 
dotted line. 

It is possible to modify the thermal and mechanical 
processing of 2024 to achieve behavioral changes in 
the desired direction, although slight alteration of 
the composition may be necessary to obtain optimum 
combinations of yield strength and tear resistance. 
The metallurgical factors involved are based on 
selective precipitation. Such changes are possible in 
any of the metallic systems that employ controlled 
precipitation for strength gains. Fortunately, con- 
siderable basic information is now available about the 
precipitation mechanisms. But the determination of 
absolute changes in yielding and tear resistance 
depend on carefully planned experiments. The near 
future will see some results from present efforts in 
this direction, as commercial advantages from suc- 
cessful use of these concepts can be most rewarding. 
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Lunar Tidal Variations in the lonospheric Layers* 


N. NARAYANA RAO 
Graduate Student in Electrical Engineering 


It has been known for 
many years that the iono- 
spheric layers are subject to 
tidal variations. In order to 
explain how lunar tidal vari- 
ations affect the ionosphere, 
we need to know the mean- 
ing of the word tide and 
to understand the behavior 
of the ionosphere. Although 
tides are commonly associ- 
ated with the oscillatory mo- 
tion produced in the ocean 
by the gravitational pulls of 
the moon and the sun, the same term is applicable 
to the atmospheric movements. However, there is a 
fundamental difference between ocean and atmos- 
pheric tides. The movements of the water are mainly 
horizontal, whereas in the atmosphere the movements 
are more complicated, due to the compressibility of 
the gases. The tidal motions can have different veloci- 
ties, both horizontal and vertical, and even different 
directions at different heights. The effect of the 
earth’s magnetic field further complicates the tidal 
phenomenon in the ionosphere. 

The term ionosphere applies to the region of the 
upper atmosphere extending from a height of about 
60 km up to about 500 km above the earth. As the 
name implies, the region consists of ionized gases, 
primarily caused by the solar ultraviolet radiation. 
The density of ionization is not uniform with height, 
and there are distinct layers with corresponding 
maximum ionizations. These are named D, E, F,, and 
F, layers. Figure 1 shows the approximate locations 
of these layers and the corresponding electron densi- 
ties. Besides these normally regular layers, a spor- 
adic E layer (known as E,) frequently exists. These 
ionized layers often behave like reflectors of radio- 
frequency energy, by bending the waves and return- 
ing them to earth through a process of refraction. 
However, if the frequency of the wave is increased 
beyond a certain limit, depending upon the value of 


N. N. Rao 


* This article is based on Mr. Rao’s thesis submitted as 
partial fulfillment of the requirements for the degree of M.S. 
in Electrical Engineering. The authors wish to express their 
thanks to the Engineering Experiment Station of the Uni- 
versity of Washington for the financial assistance rendered 
page Ss course of the investigation upon which the article 
is based. 
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maximum electron density 
of a layer and the angle of 
incidence of the wave upon 
the layer, the wave may 
penetrate through the layer. 
The highest vertically inci- 
dent wave frequency that 
can be reflected from the 
height of maximum ioniza- 
tion density of any layer is 
defined as the “critical fre- 
quency” of that layer. The 
critical frequency is hence 
dependent upon the value 
of maximum ionization density of the layer. 

The purpose of this paper is to evaluate theoretic- 
ally the effects of lunar tides upon the ionospheric 
parameters and to compare the theoretical with the 
experimental results. 


H. M. Swarm 


IONOSPHERIC DATA 


The following nomenclature is appropriate for the 
ionospheric parameters, some of which are obtainable 
by using a sweep-frequency ionospheric recorder of 
the type described in the July, 1958, issue of Trend 
by Swarm et al’: 


h'E(E,, F, and F,) =the minimum equivalent 
height of reflection of the 
E(E,, Fi, and F;) layer, 


=the height of maximum 
ionization of the 
E(F; and F,) layer, and 


=the critical frequency of the 
E(F; and layer. 


In brief, the method of collecting data is to send 
a pulse of radio-frequency energy vertically upwards 
and measure the time delay between the transmitted 
pulse and the pulse received after reflection from the 
ionosphere. If this time delay is multiplied by the 
velocity of light, the result gives the “equivalent” or 
“virtual” height at which the pulse is reflected back 
to the earth. This is not equal to the “actual” or true 
height, since the velocity of propagation of the pulse 
through the ionized region is reduced, especially as 
the incident wave frequency approaches the critical 
frequency of the layer. By varying the frequency of 
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the transmitted pulse, a virtual height-vs-frequency 
record—an ionogram—may be obtained. Figure 2 is 
a representation of such an ionogram, redrawn for 
clarity. As is apparent from the figure, the minimum 
equivalent heights of reflection and the critical fre- 
quencies can be directly ascertained from the iono- 
gram. The approximate true heights of maximum 
ionization, hmax, are equal to the virtual heights read 
off the ionogram at a frequency 0.834 times f,, ac- 
cording to Booker and Seaton’. 


THE TIDAL PHENOMENON 


The distribution of ionization of the ionospheric 
layers undergoes regular diurnal, seasonal, and sun- 
spot cycle variations, and these are quite different at 
different geographical locations. All these effects are 
largely a result of changes in the ionizing radiation ; 
but to a lesser extent they are also due to tides. In 
fact, some of the peculiarities of the F, layer are 
explainable by the tidal phenomenon, and this is one 
of the main reasons for intensive research in this 
field. The tidal effects are concealed in the variations 
of h’, hmax, and f, of the ionospheric layers, which 
are the sum of contributions owing to the sun and 
moon. Fortunately, the lunar tidal variations—unlike 
the solar ones—are easily separable. The reason is 
that the periods of lunar variations average 24 hr, 
51 min and its simple submultiples, whereas the per- 
iods of solar variations coincide exactly with those 
already present in the ionospheric data, from the 
simplest considerations of electron production. The 
semidiurnal lunar variation of period 12 hr, 25.5 min 
has been observed to possess the strongest amplitude. 

The lunar tidal phenomenon in the ionosphere may 
be obtained quantitatively by calculating the ampli- 
tudes and phases of semidiurnal sinusoidal variations 
(of period 12 hr, 25.5 min) in the ionospheric para- 


16 


VIRTUAL HEIGHT-KM 


fE 


Fic. 2. A TYPICAL IONOGRAM 


meters h’, hmax, and f,, after the solar effects have 
been removed. This process has been described by 
Swarm e? al’. Several research workers have applied 
this method to the experimental ionospheric data at 
different geographical locations. The results of these 
studies, which mostly deal with h’E, h’E,, hmaxF:2, 
and f,F, variations, are shown in Tables I and II. 
The nomenclature is as follows: 


P2 =amplitude of the semidiurnal lunar varia- 
tion in the height, 


te =phase of the semidiurnal lunar variations 
(after lunar transit), and 


200 P, =amplitude of the semidiurnal lunar varia- 
f.F, tions in the maximum electron density of 
the F, layer. 


For convenience, the critical frequency variations 
have been presented as percentage variation in the 
maximum electron density. 

So far, the presentation has dealt with the evi- 
dence of lunar tides in the ionosphere from an 
experimental researcher’s point of view. If we now 
consider the lunar semidiurnal tidal theory and its 
consequences on the ionospheric parameters, we can 
compare the results of the theoretical treatment with 
the experimental results, and can evaluate any dis- 
crepancies between the two. 


IONOSPHERIC TIDAL THEORY 

The foundation for the ionospheric tidal theory 
was laid in 1882, even before the experimental dis- 
covery of the ionosphere, when Balfour Stewart* 
suggested that certain terrestrial magnetic variations 
were an indirect manifestation of the tidal motions 
in the upper atmosphere. Stewart’s theory has become 
widely known as the “dynamo theory” of magnetic 
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TABLE I 


HARMONIC COEFFICIENTS OF SEMIDIURNAL 
LUNAR VARIATIONS IN F.-REGION DATA 


o F; 2 
GEOMAGNETIC 
LocaTION RESEARCHER MontTHS 200 MontTHS 
LATITUDE P.(km) h(br.) | of Data te(hr) | of Data 
Huancayo 1°S Martyn FF 8.4 30 3.4 4.3 36 
Ibadan 10° N Brown 7.7 8.1 40 wie 4.3 40 
Calcutta 12°N Baral 2.0 4.0 36 
Kihei 21°N Martyn 0.1 14 
Cape York 21°S Martyn x 6.4 6 af 9.4 36 
Durbanville 33° S Martyn 2.6 11.8 8 
Brisbane 36°S Martyn 3.6 $5 27 2.8 9.3 30 
Washington 41°S Martyn 2.0 6.9 84 3.3 10.4 112 
Canberra 44°S Martyn 3.7 Pe 36 a7 9.4 36 
Christchurch 48°S Martyn Bee ee 2.6 8.5 28 
Washington 50° N Martyn 11.0 12 
Hobart 52°S Martyn 2.2 7.8 9 2.4 9.3 10 
Slough 55° N Appleton and 2.0 6.0 esis aa 11.0 12 
Benyon 
Ottawa 56° N Martyn 11.4 30 
TABLE II 
HARMONIC COEFFICIENTS OF SEMIDIURNAL LUNAR 
VARIATIONS IN E-REGION DATA FOR VARIOUS SITES 
LocaTION = INVESTIGATOR LAYER P:(km) te(hr) 
Huancayo 1°S Matsushita E, 13 0.19 ice 
Ibadan 10° N Wright and Skinner E, 12 0.22 7.0 
Brisbane 36° S Martyn E and E, ay 0.50 4.5 
Brisbane 36° S Matsushita E, 13 0.70 5.7 
Brisbane 36° S Thomas and Swenson E, 10 0.69 6.7 
Brisbane 36°S Gazzard Exs 22 1.50 1.5 
Canberra 44°S Martyn E 0.19 
Stanford 44°N Swarm and Helliwell E, 6 0.59 6.4 
Washington 50° N Matsushita E, 13 0.46 5.8 
Slough 55° N Appleton and Weekes E 11 0.93 11.2 
Ottawa 56° N Datars E sa 0.40 9.5 
Ottawa 56° N Jelly E 12 0.50 0.5-5.9 


variations. According to this, the daily solar and 
lunar magnetic variations are caused by horizontal 
tidal movements of the conducting upper atmosphere 
across the magnetic field of the earth. The conduct- 
ing ionized region of the upper atmosphere may be 
pictured as the armature of the atmospheric dynamo, 
the earth being its permanent field magnet and the 
atmosphere carrying the conducting layer, the rotor. 
If a conductor moves with a velocity v in a direction 
at right angles to a magnetic field of intensity B, then 
an emf E is induced in the conductor, given by 
E = Bv per unit length of the conductor, where E, 
v, and B are components along the positive direc- 
tions of x, y, and z axes of a right-handed coordinate 
system. The induced electromotive forces, perpen- 
dicular to the earth’s magnetic field, then direct the 
positively charged particles in the opposite direction 
to the negative particles and set up a system of elec- 
tric currents, which in turn are responsible for the 
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terrestrial magnetic variations. 

Early recognition was given to the dynamo theory 
by Schuster* and Chapman’ in their work on the ter- 
restrial solar and lunar diurnal magnetic variations. 
However, when the ionosphere was discovered ex- 
perimentally, about 1925, by Appleton and by Breit 
and Tuve, the theory was confronted with serious 
quantitative discrepancies, in the conductivity of the 
ionosphere and in the velocity of the horizontal tidal 
winds. These difficulties were, to some extent, more 
apparent than real, since in calculating the electro- 
motive forces caused by tidal motions the wind veloc- 
ities have been deduced from the observed ampli- 
tudes of the barometric oscillations near the ground. 
In 1936, when Taylor® and Pekeris’ showed that in 
the general case of temperature gradient the wind 
velocities in the upper atmosphere could be much 
larger than those at the ground, the difficulties were 
only partly removed. 
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In 1932 Cowling* pointed out that intense polari- 
zation fields can be established by the currents, in 
suitable circumstances. These fields can play a more 
important part than the original induced electromo- 
tive forces in determining the ultimate pattern and 
strength of the current flow. Then it became appar- 
ent that the polarization fields can lead to a substan- 
tial increase in the effective conductivity, which re- 
lates the ultimate currents to the primary emf. On 
this basis, Stewart’s dynamo theory was subject to 
modifications by Baker and Martyn®?®", and inde- 
pendently by Hirono’? and Fejer’*®. The remainder of 
this paper extends the work of Baker and Mar- 
tyn®**"!, Their modification of the classical dynamo 
theory will henceforth be referred to as the “revised 
dynamo theory.” 

The revised dynamo problem assumes the density 
of ionization in the upper atmosphere to be inde- 
pendent of geographical position, the density being 
a function of height only, and thus treats the iono- 
sphere as a series of spherical shells surrounding the 
earth. The anisotropic conductivity of the ionosphere 
with regard to height is included in the theory. The 
geographical latitude variation is simplified by divid- 
ing the region between the geomagnetic equator and 
the geomagnetic poles into two zones of isotropic 
conductivity. One of these zones is an equatorial 
strip 3° wide on either side of the geomagnetic equa- 
tor that has been found to be effectively a better con- 
ductor than the other zone. The increased conductiv- 
ity of the equatorial zone is effective in the east-west 
direction, as has been confirmed'* experimentally. 
The height variation of the conductivity, as calculated 
by Baker and Martyn, indicates the ionosphere to 
be a better conductor at E-region heights than at 
other heights. The height of maximum conductivity 
was found to lie between 100 km to 110 km above 
the earth’s surface. Results of experimental observa- 
tions with rocket-borne magnetometers obtained by 
Singer, Maple, and Bowen ** ** and by Cahill, Jr.,’* 
indicate the existence of electric currents at 100-110 
km above the geomagnetic equator. Thus it is clearly 
established that the E-region of the ionosphere is the 
main conducting region or the dynamo layer. 


FIELDS AND CURRENTS IN THE IONOSPHERE 


The electric fields and electric currents in the dy- 
namo layer resulting from lunar tides are derivable 
by the application of the revised dynamo theory to 
a semidiurnal tidal wind-velocity potential that is 
independent of height and has the form 


sin’6 sin 2¢, (1) 


where & is an amplitude coefficient, a the radius of 
the dynamo layer, @ the geomagnetic colatitude, and 
¢@ an abbreviation for t + A + a, in which ¢ is the 
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Fic. 3. CURRENT SYSTEMS AT THE EQUINOXES THAT 
COULD PRODUCE NEW-MOON DAILY VARIATIONS 


Greenwich time in angular measure, A the longitude, 
and a a phase constant. However, the electric fields 
and currents obtainable in this way could only be 
expressed in relative magnitudes, since the value of 
k has never been estimated theoretically. This con- 
stant can be evaluated by making use of the upper- 
atmospheric current systems deduced by Chapman 
and Bartels'’ from the daily magnetic variations. A 
brief digression about the representation of these 
current systems is instructive. 

The current systems of Chapman and Bartels may 
be regarded as a world view of the relative intensi- 
ties of the daily magnetic variations at different sea- 
sons and latitudes and, in the case of lunar magnetic 
variations, at different phases of the moon. Figure 3 
depicts the current system responsible for the lunar 
daily magnetic variation on a new-moon day during 
the solar equinoxes. The currents flow in four closed 
circuits in the sunlit hemisphere, the intensity of each 
being about 5300 amperes. In the night hemisphere 
they flow in two circuits, one to the north and the 
other to the south of the equator, the intensity of each 
being about 1200 amperes. This difference between 
night and day currents conforms with the observed 
larger lunar magnetic variations during daylight 
hours than during the night. Figure 4 represents the 
mean for a lunar month, during the solar equinoxes. 
Since during a whole lunar month the daylight hours 
are equally distributed over all parts of the lunar day, 
this current system consists of eight closed circuits, 
symmetrically spaced, the current in each being about 
2600 amperes. 

Even though the Chapman and Bartels’ current 
systems were deduced from the daily lunar magnetic 
variations, it was found later that they do not fit the 
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enhanced magnetic variations near the geomagnetic 
equator. The revised dynamo theory, however, takes 
care of these equatorial enhancements. A comparison 
of the nonequatorial latitudinal variation of the Chap- 
man and Bartels’ current systems with those derived 
by the application of the revised dynamo theory leads 
to the evaluation of the unknown amplitude coeffi- 
cient, k. Thus the absolute magnitudes of the elec- 
tric fields in the dynamo layer are those represented 
in Figs. 5 and 6. Positive phase angles in these plots 
correspond to leads, from the lunar transit, at the 
rate of 30° per one hour. As can be seen, the south- 
ward electric field reverses its phase at about 35° 
latitude. 

As mentioned earlier, these electric fields in the 
dynamo layer are a sum of the dynamo or induced 
fields and the polarization fields. The polarization 
fields are coupled from the dynamo region to other 
regions along the equipotential geomagnetic field 
lines. Thus, although these fields are generated in 
the dynamo layer, they affect other regions as well. 
To estimate the effect of the polarization fields on 
the other layers, one should know what portion of 
the total electric field in the dynamo layer is due to 
polarization. Conductivity calculations and other 
relevent considerations indicate that in the equa- 
torial strip the entire field is due to polarization, 
whereas in the rest of the region about 80% of the 
total electric field is due to polarization. 


VERTICAL IONIC DRIFT 


In the ionosphere, which is considered as a series of 
thin spherical conducting shells, the magnetic field 
is inclined to the surface of the shells. The electric 
fields and currents, which are essentially in the plane 
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of the ionosphere, cause a drift of neutral ionization. 
This drift is perpendicular to both the direction of 
current and the direction of the magnetic field, and 
necessarily has both horizontal and vertical compon- 
ents. This study, however, is concerned only with 
the vertical drift velocities, v, which are related to 
the electric fields by a formula of the type 


E, 


where E£, and E, are the southward and eastward 
electric fields in the region under investigation, r, 
and r, are quantities into which the height distribu- 
tion of electron density and conductivity are lumped, 
and H is the intensity of the earth’s magnetic field. 
The results of such a calculation are plotted in Figs. 
7 and 8. 

Figure 7 indicates the amplitudes of the vertical 
ionization-drift velocities, due to the lunar semidiur- 
nal tide, at various heights ranging from 90 km to 
300 km. Figure 8 represents the phases, in hours 
after lunar transit, of these velocities. As is at once 
apparent, at E-region heights the amplitude of drift 
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velocity is large in the vicinity of the magnetic equa- 
tor and small beyond about 30° of geomagnetic lati- 
tude. A gradual change in the phase of the velocity, 
from about 8 hours after lunar transit near the equa- 
tor to about 14 hours (equivalent to 2 hours) after 
lunar transit near the poles, occurs around geomag- 
netic latitude 35°. As can also be seen, the vertical 
ionization-drift velocity has a large height gradient 
at these heights. However, at F,-region heights the 
latitudinal variation of drift-velocity amplitudes is 
not so marked. The change in phase of v is only about 
1.5 hours, from 7 hours after lunar transit near the 
geomagnetic equator to about 5.5 hours after lunar 
transit near the poles. The height gradient of v is 
also negligible in this region. 

It follows immediately that the next step in this 
treatment is to investigate the consequences of this 
vertical ionic drift on the F, and the E ionospheric 
layers. In this step the other processes in the iono- 
sphere, such as electron production, electron loss, and 
diffusion, have also to be taken into account, and 
the combined effects of these and the tidal drift must 
be determined. These are necessary considerations 
because the experimental results of the lunar vari- 
ations contain the effects of all these processes. A 
comparison of experimental results with these theo- 
retical results is then appropriate. 
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F,-LAYER VARIATIONS 


A full-fledged investigation without an approxi- 
mation to the height distribution of the ionospheric 
processes just mentioned is very hard to accomplish. 
Hence the following approximate expression for the 
velocity, W’, of vertical movement of h,,F,, the height 
of maximum ionization of the F, layer, will be 
employed : 


W = —ae~*+be~* —ce*+. (3) 


Into this expression are lumped the effects due to 
electron production, electron loss, diffusion, and tidal 
drift, which are represented by a, b, c, and v. The 
factors e** take care of the height variations. Since 
the height gradient of v is negligible at F,-region 
heights, no such factor is attached to this term. Since 
the electron production is due to solar radiation, the 
factor a becomes zero during night, whereas during 
daytime it dominates. The diffusion term has negli- 
gible influence during daytime. As a result, the effect 
of the tidal drift has to be calculated separately for 
day and night, taking into account the light-to-dark- 
ness inequality of the current systems and hence the 
drift velocities. In both the cases, the drift velocity 
produces only these perturbations : 


Ah, F,, in the height of maximum ionization density, 
and 
AN,,F,, in the value of maximum electron density. 


The method of estimating Ah,,F, consists of first 
setting v = o in the expression for W and calcu- 
lating the equilibrium height of the layer, and then 
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determining the new equilibrium height with a value 
assigned for v. The difference in the two equilibrium 
heights will then be equal to the amplitude of Ah,,F,. 
The phase lag of Ah,,F, over the drift velocity is equal 
to half of the amplitude of Ah,,F, divided by the 
amplitude of the drift velocity. The perturbation in 
the maximum electron density, expressed as fraction 
of the maximum electron density, can then be evalu- 
ated from the expression 


NaF: 
and it lags the perturbation in h,,F,, by an angle 9, 
given by 


o=arc (5) 


where £ is a coefficient representing electron decay, 
» is the angular velocity of the earth’s rotation, and 
H,, the scale height of the F, layer. 

Assuming the most reliable values ** ** ?° avail- 
able thus far for the F,-region parameters, the per- 
turbations Ah,,F, and 
function of geomagnetic colatitude for both day and 
night conditions. The solid curves in Fig. 9 repre- 
sent these calculated values of amplitudes of Ah,,F, 


are calculated as a 


and Se, to different scales. The latter is ex- 
N. 2 
pressed as a percentage so as to compare it with 


00P2 
2 
in Fig. 10 represent the calculated phase values of 
Ah,,F,. The perturbation in the maximum electron 
density lags Ah,,F, by constant amounts of two hours 


the experimentally available - The solid curves 
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during the day and three hours during the night. To 
compare these theoretical results, the corresponding 
experimental results obtained from Table I are plot- 
ted in the respective figures. From these, it can be 
easily inferred that there is, in general, a good agree- 
ment between the theoretically obtained results and 
the experimentally deduced values. 


E-LAYER VARIATIONS 

Now it will be interesting to investigate whether 
the same general agreement exists for the E-layer 
variations. For this case, the only experimentally 
available results are the lunar semidiurnal variations 
in h’E (and E,), both in amplitude and phase. Hence, 
for comparison, the effect of vertical ionic-drift vel- 
ocity on the minimum equivalent height of the E 
layer has to be evaluated. However, this evaluation, 
which requires a computer, is probably not worth 
while at this time since many of the experimental 
results are questionable, as will be indicated later. On 
the other hand, the perturbation, Ah,,£, in the height 
of maximum ionization density of E layer is easier 
to calculate by using an approximate expression, 
derived by Appleton and Lyon”', 

v 


where a is a coefficient representing electron decay 
and N,, is the maximum electron density. As can be 
seen from this expression, the height gradient of 
drift velocity does not affect Ah,,E. It should also be 
pointed out that diffusion has negligible importance in 
the E region and so was not considered in deriving 
the above expression. 
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After the values of Mitra*® and Appleton” for the 
E-region parameters are chosen, the amplitudes of 
Ah,,E for both day and night are computed at selected 
geomagnetic latitudes. The phases of Ah,,E are ob- 
tained in a manner similar to those of Ah,,F,. The 
results are presented in Table III. A comparison of 
these with the corresponding experimentally obtained 
results of h’E (and E,) variations, given in Table II, 
clearly indicates a discrepancy between the theory 
and experiment, in both amplitude and phase. 

The calculated amplitude variations from theory 
are high at low latitudes and very low at high lati- 
tudes, whereas the experimental results show an 
exactly opposite trend. The phase reversal of the 
lunar variations at about geomagnetic latitude 35° is 
not evident in the experimental results, perhaps be- 
cause most of these are for locations above geomag- 
netic latitude 35°. Even then, the phases seem to be 
scattered around 6 to 10 hours, whereas according to 
theory the maximum lunar variations at these lati- 
tudes should occur at about 0.5 to 2 hours after lunar 
transit. These discrepancies may be attributed to the 
following reasons : 

1. The published values of the minimum equivalent 
height, h’E (and E,) of the E (and E,) layer, from 
which the experimental results are deduced, have a 
very limited accuracy. Moreover, the experimental 
work at different locations was done at widely separ- 
ated intervals of time, so that the results lack con- 
sistency. This is more apparent because the results 
for the same location, but at different periods, show 
considerable divergence. 

2. For more accurate comparison between theory 
and experiment, the experimental results must be 
deduced from an analysis of actual heights rather than 
virtual heights. None of the available experimental 
results are obtained in this manner. 

3. In the theoretical treatment, the semidiurnal 
lunar-tidal wind has been assumed to be uniform with 
height in the dynamo layer. But actually it may have 
some height gradient, which cannot possibly be esti- 
mated from the presently available information. 
However, this can only slightly change the ampli- 
tudes of the variations derived from theory, and it 
does not affect their latitudinal distribution. 

4. Most of the available E-region experimental 
results shown in Table II are obtained from the 
sporadic E-layer data. According to one version of 
the causes of sporadic E ionization, irregularities in 
electron density distribution are formed in this region. 
Local polarization fields may be established at these 
irregularities; and these polarization fields might 
change the total electric fields derivable from the re- 
vised dynamo theory, and hence might contribute to 
the disagreement between the theoretical and experi- 
mental results. 
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TABLE III 
CALCULATED VALUES OF Ab,,E 
AT DIFFERENT GEOMAGNETIC LATITUDES 


GEo- DAY NIGHT 
MAGNETIC 

LATITUDE AMPLI- PHASE AMPLI- PHASE 
TUDE (km) (hr) TUDE (km) (hr) 

hg 0.17 8.00 0.60 8.70 
10° 1.00 9.30 3.45 10.00 
36° 0.07 0.25 0.21 0.95 
44° 0.03 1.45 0.11 2.15 
56° 0.07 2.10 0.23 2.80 


In conclusion, we can say that the calculated semi- 
diurnal lunar variations in the height of maximum 
ionization density and in the maximum ionization 
density of the F, region agree well with those derived 
from the experimental data for the ionospheric, 
region. But more research is needed before we can 
resolve the discrepancy between the calculated semi- 
diurnal lunar variations in the height of maximum 
ionization density of the E region and those varia- 
tions in h’E (and E,) derived from the experimental 
data for the ionospheric E (and E,) layer. 
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THE TREND IN ENGINEERING 


The Implications of Isentropic Nozzle Flow 
for a Perfect Gas 


CHARLES P. COSTELLO 
Assistant Professor of Mechanical Engineering 


The authors have long 
been concerned with the fact 
that a basic relation of clas- 
sical thermodynamics has 
always been used in an 
application that apparently 
does not satisfy the require- 
ments imposed on the basic 
relation. The problem in 
question concerns flow in 
a nozzle at constant en- 
tropy (isentropic flow), a 
problem that appears in a 
tremendous number of ther- 
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modynamic analyses. 

To determine exact relations among the thermo- 
dynamic properties of the flowing fluid, we made an 
analysis of frictionless adiabatic flow of a Perfect Gas 
in a nozzle. We found that the relations commonly 
employed in classical thermodynamics are valid, 
although the relations are based on an equation not 
strictly applicable in the case analyzed. Our investi- 
gation also shows that there is a single nozzle design 
to provide true frictionless, adiabatic flow rather than 
a family of nozzles as is implied by standard works on 
the subject. In this article we present a straightfor- 
ward approach for obtaining the relations among 
properties that must hold for flow of fluids, even for 
the ones that cannot be considered Perfect Gases. 


Nature of the Problem 


The entropy function of a system, as defined in 
classical thermodynamics, is 


as=( T (1) 


where s = entropy, Q = heat transfer to the system 
in question, T = absolute temperature at which heat 
transfer takes place, and Reversible signifies that the 
system in question must be in complete equilibrium 
with its surroundings and within itself. 

A nozzle, a device used to transform energy avail- 
able in other forms to kinetic energy, is widely used. 
In nozzles the heat transfer to or from the surround- 
ings is usually an extremely small quantity. Since 
reversible flow implies ideal frictionless flow and 
also implies that all possible energy available in other 
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forms is transformed to 
kinetic energy, reversible 
flows are highly desirable. 
As an ideal for nozzle de- 
sign, then, a reversible and 
adiabatic model is commonly 
considered, adiabatic signi- 
fying no heat transfer. The 
stipulations of adiabatic re- 
versible behavior require 
that the entropy be a con- 
stant ds = O, by the defini- 
tion of the entropy function 
just given in Eq. (1). 

In classical thermodynamics the nozzle is analyzed 
on the basis of isentropic behavior. In order to solve 
for the kinetic energy at the nozzle outlet in terms 
of pressure drop across the nozzle, a requirement is 
to have some relation among the properties of the 
fluid flowing through the nozzle. The problem of 
evaluating exit kinetic energy becomes very elemen- 
tary once some relation among properties is derived 
that holds throughout the length of the nozzle. 

It seems possible to arrive at just such a relation if 
constant entropy flow is assumed and the flowing fluid 
is treated as a Perfect Gas. The Gibbs’ relation from 
classical thermodynamics is, in part, 
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where T = absolute temperature, u = internal ther- 
mal energy, P = pressure, and v = specific volume. 

This relation holds for a unit mass of a pure fluid 
substance not subjected to external fields. For the 
equation to hold, there must be no appreciable inter- 
nal movement of a portion of the mass relative to the 
center of gravity of the mass, although the center of 
gravity itself may move. If such relative movement 
is present or if external fields are influential, terms 
must be added to Eq. (2) to account for these factors. 
In the usual analysis of nozzle flow, it is assumed at 
the outset that such factors are inconsequential, and 
Eq. (2) is employed. 

For a Perfect Gas, one that obeys the equation 


= 


Py (3) 
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NOMENCLATURE 


A =cross-sectional area for flow, ft? 

C =constant, defined by Eq. 8, ft—°F 

Cp =specific heat at constant pressure, ft-lb;/Ib °F 
Cy =specific heat at constant volume, ft-lb;/lb °F 


ge =constant of proportionality in Newton's Second Law, 
Ib/Ibe ft/sec? 


=enthalpy of fluid, ft-lb;/Ib 

=thermal conductivity, lb;/hr °F 
=ratio, Cp/Cy, dimensionless 

= mass rate of flow, Ib/sec 

=molecular mass of fluid, Ib/Ib mol 

= pressure, Ib;/ft? 

= heat transfer, ft-lby 

= universal gas constant, ft-lb;/Ib mol °R 


=signifies degrees Rankine, absolute temperature 
=entropy function, ft-lb;/Ib °R 

= temperature, °R 

= internal thermal energy of fluid, ft-Ib;/Ib 
=velocity of flow in X direction, ft/sec 

= directional coordinate, ft 

=specific volume, ft*/lb 


where R = universal gas constant and m = mole- 
cular mass of the substance, it may be shown that 


=CAT. (4) 


The specific heat at constant volume, C,, is a constant 
for the Perfect Gas. Combining (2), (3), and (4) 
with the assumption that the flow is at constant 
entropy, we obtain for the fluid as a system 


¢dT , R dy 
(5) 
This is an exact differential equation that may be 
integrated between any two states achieved by the 
flowing fluid. This provides the necessary relation 
among the fluid properties anywhere along the nozzle. 
The result of the integration, when combined with 


(3), is 


= (constant) anywhere along a nozzle, (6) 


where bat, a constant for the Perfect Gas, and 


Cp = specific heat at constant pressure, a constant 
for the ideal gas. Equation (6) is invariably used in 
the preliminary design of nozzles, to arrive eventually 
at a value for the exit kinetic energy of the fluid in 
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terms of the pressure drop from entrance to exit of 
the nozzle. 

The nature of the dilemma we now face is this: As 
explained above, Eq. (2) is valid for a fluid only in 
the absence of motion within the fluid itself. Since 
Eq. (6) is derived from Eq. (2), Eq. (6) is also 
invalid if the fluid itself has any motion within it. 
Clearly, in a nozzle there is a change in the amount 
of kinetic energy stored within the fluid itself because 
individual portions of the fluid must continuously 
accommodate themselves to the converging or diverg- 
ing walls of the nozzle passage. Seemingly, then, the 
Gibbs’ Eq. (2) should be modified to account for 
motion within the system itself by adding additional 
terms. The nature of velocity variations within the 
system will now be discussed at some length, since 
their exact nature is important. 

As we must emphasize, the velocity that would 
affect Eq. (2) is not all of the velocity a stationary 
observer would see when looking at the unit mass 
of a substance flowing through a nozzle. Equation 
(2) is still valid in accelerating systems, provided 
an observer riding along with the center of gravity 
of the system observes no motion ; that is, the motion 
that has to be excluded so as to use Eq. (2) is the 
motion of one portion of the system relative to 
another and not that of the center of gravity of the 
system. But in the nozzle flow we obviously have 
the former class of motion: an observer riding on the 
center of gravity of a unit mass flowing through a 
nozzle would observe portions of the system converg- 
ing on him or diverging from him. Nor does it seem 
safe to neglect this motion as being negligible a 
priori; intuition tells us that such velocities have 
considerable magnitude in large nozzles with high 
flow rates. 

Certainly there must be motion of one portion of 
the unit mass flowing through a nozzle relative to the 
boundaries of its own system: without it the fluid 
could not accommodate to the passage walls. If 
continuity is to be maintained, the system must also 
expand in the direction of flow relative to its own 
center of gravity. The analysis to be presented will 
be primarily concerned with this second kind (expan- 
sion of the system in the flow direction), although 
both types of motion relative to the system’s own 
center of gravity are actually accounted for. The 
velocity in the x direction will be denoted by V, and 
the velocities normal to the x direction are actually 
considered as part of the internal thermal energy of 
the fluid. This is permissible within the system to be 
defined since, under the circumstances involved here, 
internal energy and transverse molecular velocities 
are equivalent. 

How much will these motions internal to the sys- 
tem alter Eq. (6)? We set out to answer this ques- 
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tion by using neither Eq. (2) nor (6) and arrived at 
remarkable results and an insight into aspects of 
isentropic nozzle flow for a Perfect Gas. In passing 
along this insight, the authors need to present their 
analysis of this intriguing problem. 


Fundamental Approach to the Problem 


The problem may be restated in this way: given a 
perfect gas, flowing steadily, reversibly, and adiabati- 
cally in a rigid insulated nozzle, what is the relation 
among properties that must hold continuously along 
the nozzle to ensure such isentropic flow ? 

The control-volume approach will be used. That is, 
fixed imaginary boundaries will be set up, separated 
by distance dx, and flow will be allowed to enter and 
leave these boundaries. Property values of the fluid 
and characteristics of the system are shown in Fig. 1. 


No Net Heat Transfer. While no heat transfer 
from the surroundings is permitted (because of the 
specification that the nozzle is insulated), fulfilling 
the requirement of adiabatic flow makes necessary 
the assurance that no heat is transferred by conduc- 
tion through the fluid itself. This is accomplished by 
writing 


dt .@T 


where K is the thermal conductivity, assured con- 
stant; A is the cross-sectional area for flow; T is the 
temperature ; and x is the length coordinate. The left 
side of Eq. (7) is heat inflow by conduction, the 
right side heat outflow by conduction. The solution 
to Eq. (7) that guarantees adiabatic behavior is 


( 
a 
where C is a constant. 


Conservation of Energy. For the system shown in 
Fig. 1, energy is conserved if the equation 


(9) 


is satisfied, because we have already removed all 
possibilities of net heat transfer. In this equation h 
is the enthalpy of the fluid, V is the absolute fluid 
velocity, and g, is the constant of proportionality in 
Newton’s Second Law. The enthalpy accounts for 
the internal thermal energy of the system and the 
flow work done in crossing the system boundaries ; V 
is the velocity observed by a stationary observer, and, 
as pointed out above, part of this velocity is due to 
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motion of the system center of gravity, which in itself 
would not preclude use of Eq. (2). But since por- 
tions of V may be due to velocities of the fluid rela- 
tive to its own center of gravity, we will not employ 
Eq. (2) for the sake of complete generality. 


Continuity Equation. Since the mass rate of flow, 
M. is to be a constant in the nozzle, continuity re- 
quires 


M= (10) 


2 
VdV= Ea (11) 
at any point in the nozzle. Combined with Eq. (9) 
this yields 


dA 


Equation (8) may be differentiated to yield 


dA _[ @T/dx 
dT /dx Je (13) 


and, by combining (12) and (13), 
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T T+¢ 
> + dP 
Pp p+dp 
+¢ 
A A+dA 
h h+dh 
ONTROL VO! 
dx 
AV 
ee By differentiation it may be established that 
= 


C being the constant of Eq. (8). 


Momentum. To satisfy the momentum equation 
for the control volume, we write: 

(Force on left of the control volume in the x direc- 
tion) = (PA + PdA), where dA accounts for the 
element of area in the x direction at the top of the 
stationary control volume, and 

(Force on right of control volume in the x direc- 
tion) = —(P + dP) (A + dA). 

The net force, AdP, (neglecting higher-order 
terms) must be equal to the change of momentum of 
the fluid in the x direction, or 


MaV 


—AdP= 


(15) 


We now employ the Perfect Gas equation of state 
(3), which when differentiated yields 


Rj 1 Tdv 
By combining (16), (8), and (10) and by differen- 
tiating the latter two equations to obtain relations for 
dA and dV respectively, we get by tedious algebra 


R Tdv 
dT /dx ] 


The constant, C, is that of Eq. (8). 

Equation (17) provides a means of obtaining a 
value for the constant, C, in terms of properties of 
the fluid. This is necessary to complete the problem 
because C also appears in Eq. (14), which was de- 
rived from energy and continuity requirements. The 
constant C must be eliminated in order to arrive at 
a relation among properties. Equation (17) yields 

dT 
aT 


C= — (18) 
dx 


Enthalpy for a Perfect Gas. To attempt to evalu- 
ate properties of the fluid in terms of conditions that 
are usually known, we must express enthalpy in terms 
of easily measured properties. The definition of en- 
thalpy is 


h=u+Py, (19) 


which is equivalent to 
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h=u+ (20) 


for a Perfect Gas. 

On a microscopic level, internal thermal energy de- 
pends on the kinetic energy of molecules. It is the 
kinetic energy of the molecules that manifests itself 
as what is commonly called temperature. If pressure, 
volume, or other properties of the system are capable 
of causing interaction among molecules, internal ther- 
mal energy will be in part determined by such prop- 
erties. In all real substances this is the case, and 
internal thermal energy is a function of temperature 
and some other variable, say pressure. In the ideal- 
ized Perfect Gas, this is not so. Any idealized sub- 
stance that obeys the Perfect Gas Law (3) ipso facto 
has no intermolecular forces; in’ such cases, internal 
thermal energy is a function solely of temperatures. 
This concept is commonly used for Perfect Gases, 
but the proof usually rests on an assumption that a 
relation like (2) holds. We wish to keep the rela- 
tions used here perfectly general. We postulate, then, 
that internal thermal energy is a function of temper- 
ature only. The equipartition of the energy theorem 
of kinetic theory* bears out this contention. 

Equation (20), then, states that enthalpy is also 
a function only of temperature for a perfect gas. The 
relation of enthalpy to temperature is 


dh=C,AT, (21) 


where C, is the specific heat at constant pressure. 


Summary of Computing Equations 

We will now substitute (18) and (21) into (14) 
and collect the important equations used to date. The 
substitution results in 


aT 
27) 


@T /dx* _@&T/dx? 
(dv)? —dTdp +7( aT /dx dvdx dT /dx vd Tdx 
dT 
(22) 
Also we have previously used 
M (continuity), (10) 
dT _ C_ (to eliminate conduction- (8) 
dx A heat transfer) 


* Introduction to Modern Physics, Richtmeyer, Kennard, and 
Lauritsen, 5th ed., (McGraw-Hill) New York, 1955, 119. 
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Py= <r (Perfect Gas), (3) 
= [ __vdA | (energy and 
geA* A (23) 
Plus Eq. 21), 
and 
dT 
Mim aT | de 
d 


(from momentum, continuity, 


Perfect Gas, and conduction). (18) 


By use of the equations above, we have ensured 
adiabatic behavior. We have also ensured that the 
flow will be frictionless and that each lamina of the 
fluid will be in equilibrium within itself and with its 
surroundings. This has been ensured by eliminating 
frictional forces in the momentum equation and by 
writing unique values of properties for the lamina of 
fluid. The latter practice requires that the fluid be in 
internal equilibrium ; and by writing only differential 
property changes we have specified that the system 
is ‘in equilibrium with its surroundings as well. 

Since the flow is reversible and adiabatic, the en- 
tropy is constant by definition (1). The question re- 
mains, what relation among properties must hold 
along the length, +, in order that all of the equations 
listed above are satisfied ? 

Obviously boundary conditions will enter this 
problem. These need only be specification of a mass 
rate of flow, pressure drop across the nozzle, and 
conditions at the nozzle entrance, conditions that are 
normally known to the designer. Then the equations, 
when satisfied, will give the nozzle area and fluid 
temperature, specific volume, and other properties, 
all as a function of +. These variations with + will be 
the ones required to assure isentropic flow. 

Upon arriving at the equations summarized above, 
the authors were prepared to eliminate the ruinous 
a priori assumption of Eq. 6. They planned to work 
out a finite-differences solution by employing a digital 
computer ; but they discovered that Eq. 6 is, in fact, 
the only possible relation that permits simultaneous 
satisfaction of all of the equations. It is remarkable 
that such a simple specification satisfies specifications 
on equations of state, continuity, momentum, heat 
transfer, and overall energy balances. 


Implications of Isentropic Flow of a Perfect Gas 


The requirements that must be satisfied for a Per- 
fect Gas to undergo isentropic flow are as follows: 


1. Tne energy owing to expansion of portions of 


APRIL, 1960 


the unit mass of a Perfect Gas relative to its own 
center of gravity must be negligible. This is absolutely 
true, since any finite effect of such expansion would 
invalidate Eq. (6), which does actually satisfy all of 
the equations that must hold for isentropic flow. 


2. For a Perfect Gas, there is a unique nozzle that 
guarantees isentropic flow for each condition of 
mass-flow rate, pressure drop across the nozzle, and 
fluid properties at the inlet. The area as a function 
of x may be obtained from the foregoing equations 
that contain A and x factors. Usually, it is implied 
that a family of nozzies exists to give isentropic flow 
with a specified pressure drop across the nozzle. 
Strictly speaking, this is not so. The unique area re- 
quirement is derived from Eq. (7), which is based 
on the requirement that conduction-heat transfer must 
be zero. This is neglected generally because the con- 
duction-heat transfer is considered negligible a priori. 
However, the equations we have listed enable the 
construction of a nozzle ensuring the absence of heat 
transfer by conduction. 


3. Finally, and perhaps most important, we can 
write a system of equations ensuring flow at con- 
stant entropy (reversible adiabatic flow) without 
resort to such relations as (2). A system of these 
equations may be solved by trial and error to obtain 
a continuous relation among properties of the fluid in 
question. In the particular case of a Perfect Gas the 
equations are satisfied by (6). For substances other 
than Perfect Gases, the approach would differ only 
slightly. If a more complex equation of state than 
Eq. (3) is used to describe a fluid, a property rela- 
tion different from (6) would result, but the proced- 
ures employed in arriving at the result would not be 
markedly different from those employed here. We can 
safely say, however, that the result would be derived 
only at the expense of considerable labor, and that 
the problem should be programmed for a trial-and- 
error solution on a digital computer. 


STRAWS IN THE WIND 


An engineering degree is no longer a passport to 
a good job? This is the opinion of the Office of Stu- 
dent Personnel of the Cornell University College of 
Engineering. Because of the intense competition for 
engineering graduates that developed a few years 
ago, many companies have been hiring liberal arts 
graduates and technicians to fill jobs previously held 
by engineers. Cornell also reports that its five-year 
curriculum is paying dividends to its graduates in 
the form of higher-than-average starting salaries. 
Civil Engineering, Vol. 30, No. 3. March. 1960. 
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Engineering Research Dividends 


Recent theses submitted for advanced degrees from the College 
of Engineering, University of Washington, will regularly appear 
in abstract form in The Trend. Complete copies of these theses are 
available for reference at the University of Washington Library. 
Address inquiries to the Director of the Engineering Experiment 
Station, University of Washington, Seattle 5. 


CRUMBLEY, D. L., “A Study of Fuel Consump- 
tion and Travel Time Characteristics of Diesel 
Powered Vehicles,” M.S. in Civil Engineering, 
1960. 


This thesis is divided into two parts: Part A, which 
is a report of the seasonal effect on fuel consumption 
of diesel-powered vehicles; and Part B, which in- 
cludes the identification and use of some of the more 
easily measured variables used to estimate travel time 
and enroute fuel consumption. 

In the section devoted to seasonal effect on fuel 
consumption, a large part of the field work proved 
to be inadequate due to the malfunctioning of the 
basic metering unit. Therefore, alternate sources of 
information were used as a basis for this portion. 
Primarily these sources were various trucking com- 
panies’ records and local climatological reports. Data 
obtained from a controlled test vehicle was used to 
make a spot check of the accuracy of the trucking 
records. 

The second portion of this paper includes the 
identification and use of rise, fall, weight, horsepower, 
and distance to estimate travel time and enroute fuel 
consumption. These empirical variables were selected 
and equated to form two families of curves, one for 
estimation of enroute fuel consumption and the other 
for estimation of travel time. 

The facts learned by this investigation could be 
applied with significance to the field of highway eco- 
nomics and for economic transport-route feasibility 
calculations for the trucking industry. 


LARSON, J. R., “Boundary Layer Transition in 
Natural Convection,” M.S. in Engineering, 1960. 


The transition of laminar to turbulent flow in a 
natural-convection boundary layer was experiment- 
ally studied to determine if the transition might be 
correlated as a function of the Grashof number or 
Grashof-Prandtl number product. To accomplish 
this, a test apparatus consisting of a vertical copper 
pipe heated internally by electrical resistance wiring 
was submerged in a test fluid of water or ethylene 
glycol. The temperature of the ambient fluid and the 
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copper-pipe surface was kept nearly constant for a 
test run. The fluids were picked so that by varying 
their temperature the Prandtl number could be varied 
from two or three for water to over 100 for ethylene 
glycol. 

The sine-shaped wave of the initial instability lead- 
ing to transition and turbulent flow was visually 
observed by means of a concentrated arc light that 
projected parallel light by the vertical edge of the 
pipe. Flash glass or a wall served as a screen for 
the image of the boundary layer. The lowest point on 
the cylinder at which instability occurred was noted 
and measured to within about 4 in. By evaluating 
the physical properties of the fluids at the average 
film temperature, the Grashof and Prandtl numbers 
were determined at transition. 

For water the initial instability occurred at a Gras- 
hof number of 5 or 6 x 10* and a Grashof-Prandtl 
number product of 2 or 3 x 10°. For ethylene glycol 
the initial instability occurred at Grashof numbers 
from 3 to 6 x 10° and Grashof-Prandtl numbers of 
1.3 to 5.7 x 10°. 

The results showed that the transition from laminar 
to turbulent flow is not a function of the Grashof 
number alone. However, no function of the Grashof 
and Prandtl numbers was found whereby the point of 
initial instability could be exactly correlated. As an 
approximation, the initial instability may be said to 
occur at a Grashof-Prandtl number product of 3 x 
10°. 


RAO, N. N., “The Semi-Diurnal Lunar Variations 
of Iorospheric Layers,” M.S. in Electrical Engin- 
eering, 1960. 


An up-to-date review of the ionospheric “dynamo 
theory” of terrestrial magnetic variations is presented. 
Experimental evidence is collected in support of the 
location of the dynamo region at a height of 100 - 110 
km. From the revised dynamo theory of Baker and 
Martyn and the ionospheric current systems deduced 
by Chapman and Bartels from the observed lunar 
daily magnetic variation L, the vertical ionization- 
drift velocities due to the semi-diurnal lunar tide have 

(Continued on page 30) 
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Notes and Comments 


ENGINEERS IN THE NEWS 


M. M. David, Associate Professor of Chemical En- 
gineering, was the co-author with H. Glaze, Jr., for- 
mer Graduate Student, and M. G. Rao, Research 
Assistant, of papers presented January 11, at the 
Golden Jubilee Symposium, Indian Institute of Sci- 
ence, Bangalore, India. Professor David and A. F. 
Kiere, former student in Chemical Engineering, co- 
authored a paper given on December 28 at the North- 
west Scientific Association in Spokane. 


R. H. Bogan, Associate Professor of Civil Engi- 
neering, has received a new grant of $15,000 from the 
Aero-Space Division of the Boeing Airplane Com- 
pany. The grant is to be used to continue the experi- 
mental investigation of the activated sludge process 
as a means of sustenance for men during space flight, 
a study described in the January Trend. 


R. N. Kusian, Director of the Environmental Re- 
search Laboratory, is leaving the University to take 
over a position in the Hazards Control Unit of the 
Lawrence Radiation Laboratory, University of Cali- 
fornia, Livermore, California. 


R. G. Hennes, Professor of Civil Engineering, has 
been assigned the task of setting up a graduate pro- 
gram in soil mechanics at the Bengal College of En- 
gineering, Calcutta, India. The work is under the 
jurisdiction of the International Cooperation Admin- 
istration. Professor Hennes will be in India for a 
year, beginning in June. 


E. M. Horwood, Associate Professor of Civil En- 
gineering, presented papers to sections of the High- 
way Research Board on January 11 and 13 in Wash- 
ington, D. C. 


J. I. Mueller, Professor of Ceramic Engineering, 
attended the annual meeting of the American Ceramic 
Society in Philadelphia during the last week in April. 
Professor Mueller is vice-president of the Society. 


James W. Souther, Associate Professor of Human- 
istic-Social Studies, has been elected Chairman of the 
English Division of ASEE. 


A University press publication, Men and Moments 
in the History of Science, edited by Herbert M. Ev- 
ans, is a collection of talks by distinguished California 
scientists delivered at meetings of a dinner club de- 
voted to the history of Science. 


The University is to have a new computer build- 
ing, which is to be completed in the summer of 1961. 
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NEW PROJECTS 


Project No. 271. Development of a Research Dem- 
onstration Flume for Open Channel Flow. Super- 
visor: E. P. Richey, Associate Professor of Civil 
Engineering ; Research Assistant: J. F. Roth, Jr. 


Project No. 272. Influence of Nuclear Radiation on 
Boiling Heat Transfer. Supervisor: P. J. Waibler, 
Associate Professor of Mechanical Engineering ; 
Research Assistant: Meldon H. Merrill. 


Project No. 273. Single-Particle Rate Studies in Ion 
Exchange. Supervisor: M. M. David, Associate 
Professor of Chemical Engineering ; Research As- 
sistant: Norman L. Eatough. 


Project No. 274. SO.-H,O-NH, Equilibria. Super- 
visors : J. L. McCarthy, Professor of Chemical En- 
gineering, and L. N. Johanson, Associate Profes- 
sor of Chemical Engineering ; Research Assistant : 
William D. Scott. 


Project No. 275. Interfacial Resistance in Liquid 
Systems. Supervisor: W. J. Heideger, Assistant 
Professor of Chemical Engineering ; Research As- 
sistant: A. Vasudev. 


Project No. 276. Reactivity Effects of Fuel Burnup 
in Nuclear Reactors. Supervisor: D. E. McFeron, 
Professor of Mechanical Engineering; Research 
Assistant : C. W. Savery. 


Project No. 277. Control-Rod Studies for UWTR. 
Supervisor : D. E. McFeron, Professor of Mechani- 
cal Engineering ; Research Assistant: J. H. Smith. 


Project No. 278. Investigation of Stresses in Piles 
and in Surrounding Soil when Piles Are Laterally 
Loaded. Supervisor: M. I. Ekse, Professor of Civil 
Engineering ; Research Assistant: J. J. Furney. 


Project No. 279. Noise Characteristics of Parametric 
Amplifiers. Supervisor: J. L. Bjorkstam, Assist- 
ant Professor of Electrical Engineering ; Research 
Assistant : R. E. Oettel. 


Project No. 280. Construction and Electrical Prop- 
erties of Tunnel Diodes. Supervisor: J. L. Bjork- 
stam, Assistant Professor of Electrical Engineer- 
ing; Research Assistant : C. Sienkiewicz. 


Project No. 281. Self-Adjusting Control Systems. 
Supervisor: R. N. Clark, Associate Professor of 
Electrical Engineering ; Research Assistant: P. C. 
Wheeler. 


Project No. 282. Afterburner-Type Smog Elimina- 
tor. Supervisor: J. C. Firey, Associate Professor 
of Mechanical Engineering; Research Assistant : 
D. Ryder. 
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Project No. 283. On the Accuracy of the Birefrin- 
gent Coating Method of Experimental Stress An- 
alysis. Supervisor: A. S. Kobayashi, Assistant 
Professor of Mechanical Engineering; Research 
Assistant: C. Larson. This project has been dis- 
continued. 


Project No. 284. Uses of a Digital Computer in En- 
gineering Design. Supervisor: P. L. Balise, As- 
sociate Professor of Mechanical Engineering ; Re- 
search Assistant: R. F. Poultney. 


Project No. 285. Measurement of Noise in a Tunnel 
Diode. Supervisor: K. Shimada, Assistant Pro- 
fessor of Electrical Engineering; Research Assist- 


ant: W. H. Eckton, Jr. 


Project No. 286. Stress Analysis of Truncated Cones. 
Supervisor: E. H. Dill, Associate Professor of 
Aeronautical Engineering ; Research Assistant: P. 
J. Chen. 


Project No. 287. Kinetic Studies of the Nitric Acid- 
Nitrous Acid-Water System. Supervisor: A. L. 
Babb, Associate Professor of Chemical Engineer- 
ing; Research Assistant: D. M. Waldorf. 


Project No. 288. A Study of the Factors Controlling 
Anaerobic Reaction Rates. Supervisor: R. H. 
Bogan, Associate Professor of Civil Engineering ; 
Research Assistant : K. K. Chin. 
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THE SEMI-DIURNAL LUNAR VARIATIONS 
(Continued from page 28) 


been calculated by the author for different ionospheric 
heights. It is found that the vertical drift velocity 
has a large height gradient at E-region heights, con- 
trary to earlier assumptions. 

The effect of this vertical drift velocity on both 
the daytime and nocturnal ionization of the E and F, 
layers is investigated. The theoretical results ob- 
tained by the author are then compared with the 
available experimental data for the harmonic coeffi- 
cients of the semi-diurnal lunar oscillations of the 
ionospheric layers. The theoretical results for the 
F,-layer variations seem to be in good agreement 
with the experimental results, whereas in the case of 
E-layer variations the comparison indicates a general 
disagreement between the two. An attempt is made 
to explain this disagreement. 


“The Use of Metallic Thermo- 
” M.S. in Mechanical 


STRALEY, R. L., 
couples as Heat Engines, 
Engineering, 1960. 


The concern of this thesis is the theoretical anal- 
ysis of the metallic thermocouple as a heat engine and 
an experimental analysis of the effect of surface-dif- 
fused impurities. 

A review of the fundamental theory of the thermo- 
electric effect and a thermodynamic analysis of the 
metallic thermocouple are preserted to determine 
their parameters. 

The variation of the thermoelectric parameters with 
the composition of binary alloys is presented. The 
general effect of a radiation field on the thermoelec- 
tric properties of a metal is discussed, and the effect 
of “thermal” neutrons upon copper is analyzed. 

An experimental analysis is made of the variation 
of the thermoelectric properties of copper resulting 
from the diffusion of nickel radially into a test rod. 
These results are compared with a theoretical anal- 
ysis. 


Why not a two-way street? The question of how big a place the humanities should occupy in engineer- 
ing education has frequently been investigated in Trend editorials; and the topic is still wide open for 
discussion. An angle of the subject examined recently in Saturday Review is sometimes neglected. In an 
article predicting the directions technology may take in this decade, Roger Burlingame says: “If we get 
time to think, we can do something with education. There is much talk about introducing studies of the 
humanities into engineering curricula and there are occasional instances of success. But how about the 
converse of this program? How about getting more science and technology into the teaching of the 
liberal arts? We should begin this at the elementary level. Surround the child with models of what he 
must live with; show him how each thing evolved from something else; show him what is behind every 
button he pushes, every picture on the screen, how 150 horses got under the hood of Daddy’s car. The 
man or woman emerging from this kind of instruction will no longer be at the mercy of a managerial 
priesthood which is coming to possess all the ‘knowhow.’ Such an educational equipment, coupled with 
natural—not synthetic—human desire, will enable us to lead the technology we now follow.” From Roger 
Burlingame, “The Sixties: Coming up in Technology,” Saturday Review, Vol. XLIII, No. 2, Jan. 9, 1960. 
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A NEUTRON ABSORPTIOMETER FOR THE ON-STREAM ANALYSIS OF BORON 
(Continued from page 3) 
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Fic. 3. CALIBRATION CURVE 


measuring annulus of 13g in. had adequate sensi- 
tivity within the range of zero to 2 g per liter of 
boron. The calibration curve of the completed unit 
is shown in Fig. 3. 

Seal Pot. The function of the seal pot is to provide 
a liquid pressure seal between the process and oper- 
ating gallery. The seal pot will always contain suffi- 
cient solution to provide 5 ft of water head in event of 
a positive pressure in the process. The process vessel 
is equipped with a pressure relief for 2 ft of water. 
Flow tests of the seal pot and sensing cell indicated 
the unit will pass 150 gal per minute with a full 
make-up tank. 


Instrumentation 


The electronic instrumentation selected for thermal 
neutron measurements is standard counting instru- 
ments that have good reliability and accuracy. A 
block diagram of the instrumentation is shown in Fig. 
4. The BF, neutron detector tube is an ionization 
chamber filled with enriched BF, gas. If a neutron 
penetrates the counting portion of the tube and is 
captured by a nucleus of the gaseous boron-10, an 
alpha particle is emitted and the boron-10 is con- 
verted to Lithium-7 : 


+ . 


The alpha particle creates an ionization, and a voltage 
pulse is produced that may be counted by suitable 
electronic circuitry. 

Immediately above the BF, tube is located a linear 
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Fic. 4. BLOCK DIAGRAM OF INSTRUMENTATION 


preamplifier. It is as near as possible to the tube in 
order to minimize loss of pulse height resulting from 
the voltage pulse from the tube being attenuated by 
the capacitance and low impedance of the cable to the 
preamplifier. The preamplifier gives a voltage gain 
of approximately 10; it has a low output impedance 
that permits counting instrumentation to be located 
several hundred feet away from the basic sensing 
unit. 

A well-stablized high-voltage power supply pro- 
vides approximately 1500 v, dc, to operate the BF, 
tube in its “plateau” region. 

A linear amplifier receives the voltage pulses from 
the preamplifier and amplifies them to a usable level. 
The pulses are then fed to the count-rate meter, which 
measures the instantaneous rate of counting and pro- 
duces a 0 to 50 mv signal proportional to the rate. 
A circular chart recorder is used to record the count- 
ing rate. An adjustable alarm circuit is located inside 
the recorder to operate the alarm lights when the 
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Fic. 5. THE NEUTRON ABSORPTIOMETER 


counting rate of the instrument exceeds a predeter- 
mined value. The time constant of the instrumenta- 
tion totals 5 sec. 


Conclusion 

A neutron absorptiometer developed for the on- 
stream analysis of boron in the range of 0 to 2 g per 
liter is shown in Fig. 5. The unit has been calibrated 
to measure boron in solutions of water, 60% nitric 
acid, 50% sodium hydroxide, and 38% sodium ni- 
trate. The monitor could be used to measure the 
concentration of any solute with a large thermal 
neutron cross section compared to the solvent. The 
construction materials must be such that they pre- 
sent a low neutron absorption cross section immed- 
iately surrounding the neutron source and monitored 
solution. Future applications might use lucite, poly- 
ethylene, aluminum, zirconium, and stainless steel as 
construction materials. 

The monitor and instrumentation have been instal- 
led in the plant and are presently being evaluated. 
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WATER DISINFECTION AT A BEACH 
(Continued from page 10) 


their apprehension about water quality. For this rea- 
son the report recommends continued operation of 
the chlorination facility until public confidence in our 
bathing beaches can be restored through the many 
sewerage improvements now being planned. 

The operation of the chlorination facility must be 
in the hands of a responsible person employed full 
time for the work. This person can perform other 
duties, but his primary duty must be this facility. If 
operation is not continuous and effective, it is a waste 
of money, and it will not engender confidence from 
the public. 

Daily records of the operation should be kept as a 
protection to the City and for the valuable data they 
provide in increasing our general knowledge of such 
installations. 

Arrangements should be made with the Seattle- 
King County Health Department for collection and 
analysis of weekly bacteriological samples from the 
West and East Beaches and from the Lake itself out- 
side the swimming areas and away from the shore- 
line. These samples can be collected by the facility 
operator. 

An automatic control should be placed in the chlo- 
rination system, a control that will shut off the chlo- 
rinator and the pump if the pump loses its suction 
and ceases to deliver water. 

For an effective bacterial kill, the chlorine residual 
in the water where bathers and swimmers congregate 
should be maintained between 0.4 and 0.6 ppm. Bac- 
teriological tests will aid in setting the desired chlo- 
rine residual. Possibly a lower chlorine residual may 
be maintained in the summer of 1960. 
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